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Abstract: We present the double real radiation corrections to the hadronic tt production 
stemming from partonic processes with fermions only. For this purpose, we extend the 
NNLO antenna subtraction formalism developed originally for the computation of jet ob- 
servables in e'^e~ annihilation to include the evaluation of hadronic observables involving a 
massive pair of particles. In all partonic processes, we checked the validity of our subtrac- 
tion terms given for leading and subleading colour contributions numerically by showing 
that the ratio between real radiation matrix elements and subtraction terms approaches 
unity in all single and double unresolved configurations. 
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1. Introduction 

The production of top-antitop pairs is one of the core measurements at LHC. Top quarks 
[1,2] are measured through their decay into a bottom quark and a subsequently decaying 
VF-boson, yielding up to six-jet final states for top quark pair production. With a mass 
mt = 173 =b 1.3 GeV, the top quark is the heaviest Standard Model particle produced 
at colliders and due to its very large mass it decays before it hadronises. By studying its 
properties in detail, it is hoped to elucidate the origin of particle masses and the mechanism 
of clcctroweak symmetry breaking. Since its discovery at the Fermilab Tevatron, a number 
of these properties (mass, couplings) have been determined to an accuracy of ten to twenty 
per cent. With the large number of top quark pairs expected to be produced at the 
LHC, the study of the top quark is becoming precision physics. This large production 
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rate aims to obtain precise measurements of its properties and production cross sections 
with an expected experimental accuracy of the order of five per cent. In order to extract 
fundamental parameters like couplings and masses by comparing theoretical predictions 
and data, it is therefore mandatory to have theoretical predictions with a comparable 
accuracy. 

In order to match the expected experimental precision for these heavy quark production 
processes at LHC, fixed-order calculations need to be considered at least at the next-to- 
leading (NLO), if not even at next-to- next-to-leading order (NNLO) in perturbative QCD. 

Most of the calculations involving top particles are performed in the so-called narrow 
width approximation. The top quark appears as virtual particle in hadron collider processes 
and due to the small ratio between its width and its mass, it is possible to factor the cross 
section of processes involving top quarks into the product of the production cross section 
for on-shell top quarks and the corresponding decay width. In this paper, we shall also 
follow this approach, in which the top quark pairs arc considered to be produced on-shell. 

Current theoretical predictions for the top quark pair production cross section include 
NLO corrections [3-5] and next-to-leading-logarithmic resummation (NLL) [6]. More re- 
cently even the NNLL resummation effects have been completed in [7]. These predictions 
lead to a theoretical uncertainty of the order of ten per cent. The same precision is available 
for single top quark production [8], top-pair-plus-jets production [9,10] and for top-pair- 
plus-bottom-pair production [11,12]. 

At present, a full fixed order calculation of the total top-antitop rate at NNLO is still 
missing. In such a calculation, three essential ingredients enter: double real dcr^^, mixed 
real- virtual, da^^ , and double virtual contributions da^^ . 

Recent progresses has been accomplished concerning the two- loop contributions. The 
parts of these two-loop virtual corrections that are built with products of one-loop virtual 
amplitudes have been already computed in [13]. Concerning the two- loop virtual correc- 
tions which are built with products of two-loop and tree-level amplitudes, the situation is 
different since the two-loop virtual corrections for the processes qq — > tt and gg tt are not 
fully available at present. However, a purely numerical evaluation of the quark-initiated 
process [14] could be partly confirmed by analytical results [15], which were most recently 
extended also to the gluon-induced subprocess. 

Most recently, progress towards the computation of mixed real-virtual corrections has 
been achieved, and the infrared structure of these contributions has been studied [16]. 
Also, following a semi-numerical approach [17] partonic processes contributing to double 
real radiation to top quark pair production have been evaluated [18]. 

In this paper, we shall evaluate the double real contributions to hadronic top quark 
pair production arising from partonic processes involving only fermions. 

Generally, for hard scattering observables, the inclusive cross section with two incoming 
hadrons Hi , H2 can be written as 



da 



a,o 



(1.1) 
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^1 and ^2 are the momentum fractions of the partons of species a and b in both incom- 
ing hadrons, fi being the corresponding parton distribution functions, ddab denotes the 
parton-level scattering cross section for incoming partons a and h which depends on the 
the renormahsation and factorisation scales denoted by jiR and fip respectively. 

The partonic cross section d&ah has a perturbative expansion in the strong coupling 
which depends on the renormahsation scale fXR. For a hadronic observable, a theoretical 
prediction is obtained at a given order in when all partonic channels contributing at 
that order to the partonic cross section are summed and convoluted with the appropriate 
parton distribution functions as in eq.(l.l). 

Beyond leading order, these partonic channels contain both ultraviolet and infrared 
(soft and collinear) singularities. The ultraviolet poles are removed by renormahsation, 
while collinear poles originating from the radiation of intial-state partons are cancelled by 
mass factorisation counterterms. The remainining soft and collinear infrared poles cancel 
among each other only when all partonic channels are summed over [19]. 

While infrared singularities from purely virtual corrections arc obtained immediately 
after integration over the loop momenta, their extraction is more involved for real emis- 
sion (or mixed real-virtual) contributions. There, the infrared singularities only become 
explicit after integrating the matrix elements over the phase space appropriate to the dif- 
ferential cross section under consideration. Since these obscrvables depend in general in 
a non trivial manner on the experimental criteria needed to define them, they can only 
be calculated numerically. The computation of hadronic observables including higher or- 
der corrections therefore requires a systematic procedure to cancel infrared singularities 
among different partonic channels before any numerical computation of the observable can 
be performed. Subtraction methods explicitly constructing infrared subtraction terms are 
well-known solutions to this problem. 

The crucial points that all subtraction terms must satisfy are that (a) they approximate 
the full real radiation matrix elements in all singular limits, (b) they are still sufficiently 
simple to be integrated analytically, (c) they should account for the limit they are aimed 
at without introducing spurious infrared singularities in other limits. 

For massless QCD, and for the task of next-to-leading order (NLO) calculations [20] 
of hadronic and jet observables, various subtraction methods have been proposed in the 
literature [21-25]. The dipole formalism of Catani and Seymour [21] and the FKS [22] of 
Frixione, Kunszt and Signer are the most widely used ones. Both have also been imple- 
mented in an automated way in [26-31]. 

Several methods for systematically constructing general subtraction terms have been 
proposed at NNLO [17,18,32-46]. Another NNLO subtraction scheme has been proposed 
in [47]. It is not a general subtraction scheme, but it can nevertheless deal with an entire 
class of processes, those without coloured final states, in hadron-hadron collisions. It has 
been explicitly applied to several observables [47-52]. 

In addition, there is the sector decomposition approach which avoids the need for 
analytical integration, and which has been developed for virtual [53-56] and real radia- 
tion [56-60] corrections to NNLO, and applied to several observables already [61-64]. 
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In order to evaluate double real contributions to tt production at hadron colliders, 
we will follow the NNLO antenna subtraction method. This method has been originally 
developed for the production of massless partons in electron-positron annihilation [32, 66, 
67], leading to the successful description of the infrared structure of three-jet events at 
NNLO [68,69] and the subsequent numerical calculation of the NNLO corrections to event 
shape distributions [70-73], the moments of event shapes [74,75] and jet rates [76,77]. 

For processes with initial-state partons and massless final states, the antenna subtrac- 
tion formalism has been so far fully worked out to NLO in [78,79]. It has been extended to 
NNLO for processes involving one initial state parton relevant for electron-proton scatter- 
ing in [80]. Several essential steps towards a full NNLO calculation of jet rates at hadron 
collider have been achieved in [81-84]. 

For QCD observables involving massive particles, subtraction methods have been so 
far only developed completely up to the NLO level using generalisations of the dipole 
formalism [85-87], or of the antenna formalism [88,89]. An extension of the antenna 
formalism to deal with the production of a massive pair of particles in association with jets 
up to NNLO will be presented in this paper. 

For these observables, QCD radiation from massive particles can lead to soft diver- 
gencies but cannot lead to strict coUinear divergencies, since the mass acts as an infrared 
regulator. As a consequence, less divergent contributions arise and the infrared structure 
of such processes is expected to be simpler. The kinematics is more involved due to the 
presence of finite values of the parton masses, though. In addition, in calculations of such 
observables, logarithms involving the ratio of the scales present in a given reaction occur. 
These finite effects are related to the singular behaviour of the matrix elements in the 
massless limit [85,86,88,89]. Depending on the kinematics of the reactions involved, these 
logarithms, although finite, may or may not be enhanced. 

For the kinematical situation under consideration in this paper, i.e. the production of 
a heavy fermion pair whose mass is of the the same order as as the partonic center-of-mass 
energy, we shall ignore these finite effects as also done in [18], since they are not expected to 
be large. As a consequence, in our calculation of the double real corrections to tt hadronic 
production, singular infrared behaviours are considered to arise only through the presence 
of collinear and/or soft radiation of massless particles and only subtraction terms capturing 
these features will be considered in the context of the present paper. 

For hadronic processes with massive fermions and jets in the final state, the antenna 
subtraction formalism has been extended in [88, 89] up to NLO. In [89] subtraction terms 
capturing single unresolved, i.e. soft or collinear, radiation of the real contributions to the 
hadronic reactions pp — > QQ and pp — > QQ + jet were explicitly constructed for leading 

and subleading colour contributions. The results obtained in this previous paper provide 
us with essential ingredients for the construction of subtraction terms for the hadronic 
reaction pp — >• QQ evaluated at NNLO and considered in this paper. 

At NNLO, any massless or massive calculation of the partonic contributions to an 
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hadronic observable involving an m particle final state is given by 

d&NNLO = / d^NNLO + / ^^N^LO + / ^^N^LO^ (l-^) 

where there are three distinct contributions due to: double real radiation da^^^Q, mixed 
real- virtual radiation dcj^^^^Q and double virtual radiation da^^j^Q, which are separetely 
infrared divergent. Employing a subtraction method, eq.(1.2) amounts to [32] 

^NNLO 



.MF,1 
^NNLO 

fm+l 



Jd^m+l Jd^m+l 

+ / da^NLO+ / da^NLO^ (1-3) 

where da^j^^Q denotes the subtraction term for the (m-|-2)-parton final state which behaves 
like the double real radiation contribution da§^i^Q in all singular limits. Likewise, daJ^'^^Q 
is the one- loop virtual subtraction term coinciding with the one- loop (m -|- l)-final state 
d^N^lLO singular limits. The two-loop correction to the (m -|- 2)-parton final state is 

denoted by da^^j^Q. In addition, as there are partons in the initial state, there arc also two 
mass factorisation contributions, da^^'^Q and da^^'^Q, for the (m -|- 1)- and m-particle 
final states respectively. 

In order to construct a numerical implementation of the NNLO cross section, the 
various contributions must be reorganised according to the number of final state particles. 



da NNLO = / [da NNLO " da'NNLo] 

+ / [d^NNLO ~ dd-NNLo\ 

+ / [d&VNLO-da^NLo]^ (1-4) 



where the terms in each of the square brackets is finite and well behaved in the infrared 
singular regions. 

Due to the presence of massive fermions in the final state, fewer subtractions terms 
are required since the real matrix elements develop singular behaviours in fewer regions of 
phase space than in the massless case. However, the kinematics is more involved due to 
the finite value of the parton masses and the integration of the subtraction terms is more 
difficult. A major step concerning the latter has recently been accomplished in [90], where 
the authors integrated one of the genuine NNLO antenna functions. We shall also use their 
results in this paper. 

More precisely, the purpose of this paper is to evaluate the first line of eq.(1.4), i.e. to 
evaluate da^^^^Q — dafjj^^Q applied to the production of tt at LHC using a generalisation 
of the antenna formalism at NNLO. Towards a computation of all partonic channels con- 
tributing to this process, we shall here restrict ourselves to the partonic processes involving 
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only initial and final state fcrmions. Those can be obtained by crossing any two massless 
final state fermions to the initial state in the fictitious process — )• QQqqq'q' , where we 
also allow the possibility that both massless quark-antiquark pairs are of the same flavour. 
It is worth noting here that the contributions coming from the identical-flavour process 
qq QQqq as well as the contributions from the process qq' QQqq' with its identical 
flavour counterpart are presented here for the first time. 

The plan of the paper is as follows: In section 2 we present a generalisation of the 
antenna subtraction formalism to evaluate hadronic reactions involving a final-state heavy 
quark pair at NNLO. There, we present the construction of the double real emission sub- 
traction terms capturing single and double unresolved contributions of real radiation matrix 
elements related to the all-fermion processes. We particularly focus on the changes intro- 
duced by the presence of massive final states in the expressions of the subtraction terms 
compared to the case in which only massless partons are present. In section 3 we give a list 
of all genuine NNLO four-parton tree-level antennae required. In section 4, after having 
presented all massless and massive double unresolved factors encountered in singular con- 
figurations of the double real matrix elements and four-parton antennae, we give the single 
and double unresolved limits of the massless and massive four-parton tree-level antennae 
encountered in our calculation. Section 5 contains our main results. There, we present 
the colour decomposition of the real contributions and give their associated subtraction 
terms for the all-fermion processes involved. In all cases, leading and subleading colour 
contributions are considered. In section 6, we test the validity of the subtraction terms. 
We check that the ratio between the real radiation matrix elements and the corresponding 
subtraction terms approaches unity in all single and double unresolved limits. Finally, 
section 7 contains our conclusions. 

Two appendices are also enclosed: appendix A gives a list of all the tree-level three- 
parton antennae required in the context of this paper, while appendix B summarises the 
phase space mappings required to construct the subtraction terms presented in section 5 
in all configurations. 

2. Antenna subtraction for double real radiation with massive final state 
fermions 

Beyond the leading order and up to the next-to-next-to leading order in perturbation 
theory, the trcc-lcvcl matrix elements squared associated to the real emission contributions, 
when integrated over the phase space, develop singularities when one or two final-state 
partons are unresolved, i.e. are soft or collinear. The construction of subtraction terms 
can face this problem as mentioned in section 1. 

The aim of this section is to present the general structure of the subtraction terms 
which capture the single and double unresolved features of the real radiation contributions 
related to the hadronic production of a pair of heavy quarks QQ in association with (m — 2) 
jets at NNLO in perturbative QCD. For this purpose, we develop an extended version of 
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the antenna formalism, in order to account for hadronic processes with a massive fermion 
pair in the final state. 

As mentioned in section 1, for the kinematical situation under consideration in this 
paper, i.e. the production of a heavy fermion pair whose mass is of the the same order as 
the partonic center-of-mass energy, the single and double unresolved features are related 
only to collincar and/or soft radiation of up to two massless partons. Mass effects, however, 
are present in the infrared limits of the matrix elements squared: The presence of massive 
partons in the final state changes the kinematics and influences the soft behaviour of the 
real matrix elements in a non-trivial way as explained in sections 4 and 5. 

Quite generally within the antenna formalism, the subtraction terms, which reproduce 
the singular behaviour of trcc-lcvcl real matrix elements, are constructed from products 
of antenna functions with reduced matrix elements. These antenna functions capture all 
unresolved radiation emitted between two hard partons, the radiators, which can be either 
massless or massive. The subtraction terms can be integrated over a phase space which 
is factorised into an antenna phase space (involving all unresolved partons and the two 
radiators) multiplied by a reduced phase space, where the momenta of the radiators and 
the unresolved radiation arc replaced by two redefined momenta. These redefined momenta 
can be in the initial or in the final state depending on where the corresponding radiator 
momenta are, and they are defined by appropriate mappings. Depending on where the 
two radiators are located, in the initial or in the final state, we distinguish three types of 
configurations: final-final, initial-final and initial-initial. They will be denoted with the 
shorthand notation: (f-f, i-f, i-i). The full subtraction term is then obtained by summing 
over all antennae required in one configuration and by summing over all configurations 
needed for the problem under consideration. 

Subtraction terms are constructed using the fundamental factorisation properties of 
QCD amplitudes in their soft and collinear limits. However, in QCD (unlike in QED) these 
fundamental properties are satisfied by colour-ordered sub-amplitudes rather than by the 
the full amplitudes. The antenna formalism uses colour ordering properties in an essential 
way: Real emission matrix elements are decomposed into a linear combination of colour- 
ordered amplitudes, and subtraction terms are constructed for each of these amplitudes 
separetely. Furthermore, in colour-ordered amplitudes, only partons which are colour- 
connected can lead to a singular behaviour of those. This helps considerably to determine 
the type of antenna functions that are to be used in each subtraction term. 

At the next-to-leading order, antenna subtraction terms are constructed solely with 
tree-level three-parton antenna functions which can be massless [32, 78, 79] or massive [88, 
89] depending on the process under consideration. These three-parton antenna functions 
encapsulate all singular limits due to the emission of one unresolved parton between two 
colour-connected hard radiators. 

At NNLO, the subtraction terms aiming to capture the single and double unresolved 
features of the double real matrix elements squared are constructed with either products of 
two three-parton antennae, or with genuine NNLO four-parton antennae. This will depend 
on the colour connection between the unresolved partons as explained in the massless case 
in [32,82] or as will be explained below in the massive case. 
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Compared to the massless antenna formalism, the presence of massive partons in the 
final states modifies the subtraction terms in a non-trivial way. Parton masses lead to 
modified kinematics and have to be taken into account for the phase space factorisations. 
Furthermore, those also modify the soft behaviour of the real matrix elements and it is 
precisely in order to properly capture these soft behaviours within this formalism, that 
antennae with massive radiators have to be considered. 

For the NNLO corrections to QQ + jets production in hadronic collisions, we will 
need all three types of subtraction terms and therefore all three types (f-f, i-f, i-i) of three 
and four-parton antenna functions, involving either massless or massive radiators. Three- 
parton antennae (massless and massive) have been derived before. Those of relevance in our 
context, the tt hadronic production from all-fermion process, will be given in appendix A. 
All required four-parton antenna functions will be presented in section 3 and their limits will 
be given in section 4. The phase space factorisation required to define the genuine NNLO 
subtraction terms, i.e. those involving four-parton-antennae will be presented below. 

2.1 Double real radiation contributions to heavy quark pair production 

We are interested in describing the double real radiation corrections to the process 

pp ^ QQ + {m - 2)iets. (2.1) 

To setup the normalisation, we first present the leading order (LO) m-partonic contribution 
to this process where the partons can either be massless or massive. It reads, 

daLo{Pl,P2) =M'lO X] ^^rn{PQ ;Pl,P2) 
m-2 

X-^^ \Mm+2{PQ,PQ,P5, ■ ■ ■ ,Pm-2;Pl,P2)\'^ J^\pQ,Pq,P5, ■ ■ ■ ,Pm-2)- (2.2) 
'-'m-2 

pi and p2 are the momenta of the initial state partons, the massive partons Q and Q 
have momenta pQ and pg, while the momenta of the remaining (m — 2) massless final 
state partons are labelled p^ . . ■Pm-2- Sm-2 is a symmetry factor for identical massless 
partons in the final state. Jrrr\pQ,PQ,P5, ■ ■ ■Pm-2) is the jet function. It ensures that 
out of (m — 2) massless partons and a pair of heavy quarks Q and Q present in the final 
state at parton level, an observable with a pair of heavy quark jets in association with 
{m — 2) jets is built. At this order each massless or massive parton forms a jet on its 
own. The normalization factor Mlo includes all QCD-independent factors as well as the 
dependence on the renormalised QCD coupling constant a^. X^^_2 denotes the sum over 
all configurations with (m — 2) massless partons. is the phase space for an m-parton 

final state containing (m — 2) massless and two massive partons with total four-momentum 
Pi +1^2- 111 d = 4 — 2e space-time dimensions, this phase space takes the form: 

d^m{pQ,PQ,P,, . . . ,Pm-2;Pl,P2) = ^^^(^^ ^E^^^' 

d*^ ^P5 d'^ ^Pm-2 fr, \d xdf , \ Co o\ 

^2^^(2^ ••• 2^^_,(2.)^-i ('"^ nPi+P2-PQ-PQ-P.-...Pm-2).i2.^) 
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In eq.(2.2) |A^m+2p denotes a colour-ordered tree-level matrix element squared with 
m final state partons, out of which two are massive and two are initial state partons. These 
terms only account for the leading colour contributions to the squared matrix elements. On 
the other hand, subleading colour contributions involve in general interferences between 
two colour-ordered parton amplitudes. However, to keep the notation simpler we denote 
these interference contributions also as |7Wm+2p- 

The NLO real radiation partonic contributions to the production of a heavy quark 
pair in association with (m — 2) jets involve (m + l)-final state partons (one more parton 
than the leading-order case) with two of them being massive. They have been derived 
in [89] together with their corresponding subtraction terms capturing all single unresolved 
radiation. 

Here we are interested in the construction of the subtraction terms for the double real 
radiation do^f^J^^lo '^^lich correctly subtract all single and double unresolved singularities 
contained in the matrix elements with (m + 2)-final state partons present in the double 
real emission contributions. These may be written as 

'^^NNLo{Pl^P2) = M'nNLO XI ^^rn+2{PQ,PQ,P5, Pm;Pl,P2) 

m+2 

X-^^ \Mm+4:(.PQ,PQ,P5, ■ ■ ■ ,Pm;Pl,P2)\'^ J^'^'^\pQ,Pq,P5, ■ ■ ■ ,Pm)- (2.4) 
'^m+2 

= J^NNLO XI d$rn-h2(P3, ■ • • , Pm+4; Pi , P2) 
m+2 

X^^\Mm+4iP3,---,Pm+4,;Pl,P2)f J^^^\p3 , ■ ■ . , Pm+4) , (2.5) 
^m+2 

where the last line is obtained by relabelling all final state partons. This procedure enables 
us to follow the masslcss notation as given in [82, 83] more closely. In this equation, the 
jet function Jm ensures that out of m massless partons and a QQ pair, an observable 
with a pair of heavy quark jets in addition to (m — 2) jets, is built. 

This contribution to the NNLO cross section develops singularities if one or two partons 
are unresolved (soft or coUinear). Depending on the colour connection between these 
unresolved partons, the following configurations must be distinghuished [32,82]: 

(a) One unresolved parton but the experimental observable selects only m jets. 

(b) Two colour-connected unresolved partons (colour-connected). 

(c) Two unresolved partons that are not colour-connected but share a common radiator 
(almost colour- unconnected) . 

(d) Two unresolved partons that are well separated from each other in the colour chain 
(colour-unconnected) . 

(e) Compensation terms for the over-subtraction of large angle soft emission 
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This separation among subtraction contributions according to colour connection is 
valid in final-final, initial-final or initial-initial configurations and, in any of them, the sub- 
traction formulae have a characteristic structure in terms of the required antenna functions. 
This antenna structure has been derived for processes involving only massless partons, for 
the final-final and initial-final cases in [32, 80] and in [82] for the initial-initial case. The 
presence of massive partons in the final state does not modify the general structure of the 
subtraction terms required to match the unresolved features given above in itemized form. 

Concerning the items related to the colour-connection of the unresolved partons, the 
configuration (c), the almost colour-unconnected case, and the configuration (e), regarding 
the treatment of large angle soft radiation, do not occur in the all-fermion processes that 
contribute to the double real corrections to ti production. On one hand, there are not 
enough final state partons (besides the pair of massive ones) to cover the configuration (c) 
and, on the other hand, the absence of final state gluons forbids configuration (e). In the 
following, we shall therefore restrict ourselves to describe how the configurations (a), (6) 
and (d) are dealt with here while leaving the discussion of configurations (c) and (e) to be 
treated elsewhere. 

2.2 Subtraction terms for single unresolved radiation da^^^Q 

The subtraction terms for single unresolved radiation associated to QQ production in 

association with (m — 2)- jets can be taken from the expressions derived in [89]. It is recalled 
below in all three final- final (f-f), initial- final (i-f) and initial- initial (i-i) configurations. 

perms '^"^+2 

X Y \Mm+3i- ..,I,K,.. .)|V(,-+1)(P3, ...,PI,PK,... ,Pm+4) , (2.6) 

j 

'^^mtio = ^NNLO Y d*m+2 (P3 , • • ■ , Pm+4 ; Pi , P2 ) 

perms '^"^+2 

E Y^^jk \MmU- ..J,K,.. .)\^J^+^\ps, . . . ,PK, . . . ,Pm+4) , (2.7) 

i=l,2 j 

~'S',a,(u) 



NNLO = ■N'nNLO Y ■^^"1+2 (P3 , • • • , Pm+4 ; Pi , P2 ) - 



Sm+2 

perms ^ 

X E E^Ll-^-+3(---'^'^'---)l'4"^"''HP3,...,Pm+4). (2.8) 
«,fc=l,2 3 

In these equations, the sum over j is the sum over all unresolved partons in a colour-ordered 
amplitude between radiators i and k which can be both located in the final state eq.(2.6), 
i in the initial state and k in the final state eq.(2.7) or both in the initial state eq.(2.8). 

In the reduced matrix elements |A^P, the redefined momenta are respectively denoted 
by (/, K) in the final-final, by (/, K) in the initial-final and by (/, K) in the initial-initial 
configurations, respectively. We shall adopt this notation throughout the paper: Initial 
state momenta will be denoted with a hat everywhere. 
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The subtraction terms given in the equations above, involve the phase space for the 
production of (m + 2) partons, with two of them being massive, the three-parton 

antenna functions appropriate to a given configuration Xfj^^, X^-^ and X^^,-, the colour- 
ordered reduced (m -I- 3)-parton amplitude squared |A4m+3p and the jet function J^'^^\ 

Compared to the NLO expressions given in [89] the NLO jet function is now 

replaced by J^"^^^ in all three configurations. It ensures that out of (m — 1) massless 
partons and a pair of massive final state partons, (m — 2) jets and a QQ jet pair is built. 
The jet function and the reduced m-parton amplitude do not depend on the individual 
momenta pj, pj and pk but will only depend on the redefined momenta pi and pK which 
are linear combinations of the original momenta Pi,Pj,Pk in the appropriate configuration. 
The three types of three-parton antenna functions Xfjj^ ,Xf -j^ and X?^ ■ can involve either 
massless or massive hard radiators. They depend only on the original momenta pi, pj and 
Pk , and in the presence of massive hard radiators, also on the masses of those final state 
partons. These antennae have all been derived before and the ones relevant in the context 
of this paper are given in appendix A for completeness. 

Depending on the configuration considered, the momenta present in the reduced matrix 
elements and in the jet function are defined with a particular 3 — >■ 2 mapping. Specially, 
for the initial-initial mapping, all momenta, massless and massive, in the arguments of the 
reduced matrix elements |A4m+3p and the jet function J^~^^^ have to be redefined. They 
are denoted with tildes in eq.(2.8). The two hard radiators are simply rescaled by factors 
xi and X2 respectively. In [89] we showed that the original initial-initial mapping derived 
in [79] for massless partons only, can also be used for processes involving massive partons. 
These mappings are recalled in appendix B. 

In the subtraction terms above, da^^^^ coincides with the matrix element (2.5) when 
j is unresolved. At NNLO, however, the jet function J^'^^^ allows one of the (m -|- 
1) momenta to become unresolved. As a consequence, the matrix element present in 
the subtraction terms can have an unresolved parton. In this limit d(T^^^Q does not 
coincide with the matrix element (2.5) and, as explained in detail in [32,82], these spurious 
double unresolved configurations will cancel against similar unresolved features encountered 
for the subtraction terms coming from other configurations. In this paper, those other 
configurations will arise only for the items (6) or (d) , associated to two unresolved particles 
which are colour-connected or colour-disconnected respectively. 

The subtraction terms <^cf^%lo defined above, need to be added back in their inte- 
grated form to the {m + 1) partonic contributions. Poles from those will cancel against 
explicit poles in the one-loop virtual do"^^ contributions and against those present in the 
mass factorisation contributions da^^'^ given in eq.(1.2). 

To define the integrated forms of these subtraction terms d(j^;",^Q involving only three- 
parton antennae, an appropriate phase space factorisation needs to be considered in all 
three configurations. Those have been derived in the massless case in [32, 79] and in the 
massive case in [88,89]. 
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2.3 Subtraction terms for colour-connected double unresolved radiation da^y^^Q 

When two unresolved partons j and k are adjacent, we build the subtraction term starting 
from the four-particle tree- level antennae in all three configurations Xijki, Xijkl and Xnjk- 
By construction, they contain all colour-connected double unresolved limits of the (m + 4)- 
parton matrix element associated with partons j and k unresolved between radiators i and 
I. However, these antennae can also become singular in single unresolved limits associated 
with j or k where they do not coincide with limits of the matrix elements. To ensure that 
these subtraction terms are only active in the double unresolved limits of the real matrix 
elements, we subtract the appropriate single unresolved limits of the four-particle tree-level 
antennae. For this purpose we use products of two tree-level three-particle antennae in the 
appropriate configurations. As in the single unresolved case, we replace the original hard 
radiators with new particles, I and L. When one of the hard radiators is in the initial 
state, pf=pi = XiPi and when both are in the initial state, pj = pi = XiPi, p^ = pi = xipi 
and all other momenta have to be Lorentz boosted. 

The colour-connected double subtraction term in all three configurations then reads: 



m+2 



perms 
jk 

X \Mm+2{...,I,L,...)\^jtH-..,Pi,PL,...), (2.9) 

^^n'nLO = ■^NNLO XI d$m+2 (P3 , ■ • • , Pm+i ; Pi , P2 ) 

perms '^'"+2 

^ X/ X/ i^tikl ~ ^i,jk^I,Kl ~ ^jkl^i,JL) 
i=l,2 jk 

X \Mm+2{...J,L,...)fjt\...,PL,...), (2.10) 

'^^n'nLO =-^NNL0 d$^+2 (P3 , • • • , Pm+4 ; Pi , P2 ) 

'-'m+2 

perms ' 

^ X X {-^il,jk ~ ^l,jk^iL,K ~ ^i,kj^Il,j) 

il=l,2 jk 

X \Mm+2i...J,L,...)\^jt\p3,...,Pm+4), (2.11) 

where the sum runs over all colour-adjacent pairs j, k and implies the appropriate selection 
of hard momenta i, I which, as usual, have three possible assignments of radiators. In all 
cases the (rn + 2)-parton matrix element is evaluated with new on-shell momenta given by a 
momentum mapping appropriate in each configuration. Those 4 — > 2 mappings have been 
derived in [79, 82] in the massless case. We checked that the same mappings, which are 
presented in appendix B for completeness, can also be applied in the presence of massive 
partons. 

The d(Tj\/^^Q subtraction terms involve: the phase space for the production of (m + 
2) partons, d<l>m+2i with two of them being massive, the colour-ordered reduced (m -|- 
2)-parton amplitude squared |A^^„+2p, the products of three-parton antennae and four- 
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parton antennae appropriate to a given configuration and the jet function Jm ■ As before, 
reduced amplitudes and jet functions only depend on redefined momenta in the appropriate 
configurations, while the tree-level four-parton antennae depend on the original momenta 
z, J, k, I. 

The relevant massless and massive four-parton antenna functions required to derive the 
subtraction terms for the double real contributions to ti production coming from fermionic 
processes will be given in section 3 while their single and double unresolved limits will be 
given in section 4. 

2.3.1 Phase space factorisation for two colour-connected unresolved particles 

As mentioned before and explained in [83] , the part of da^^^^Q involving products of three- 
parton antennae needs to be added back in its integrated form at the (m -|- l)-parton level. 

The part of da^^^j^Q involving four-parton antennae, instead, needs to be added back in 
its integrated form as an m-parton contribution. It can then be combined with the virtual 
two-loop contributions da^^ and the mass factorisation counterterms da^^'^ , which have 
also m final state partons as shown in eq.(1.3), and as explained in [32,82,83]. 

S b 

In order to be able to integrate these subtraction terms, da^^j^j^Q involving four-parton 
antennae, one therefore needs the 2 — )• {m + 2) parton phase space to factorise into an 
antenna phase space involving four-partons, the unresolved partons and the radiators, 
times a reduced 2 — t- m phase space. This factorisation will be different depending on 
whether the hard radiators arc in the initial or in the final state. For the cases involving 
only massless particles the phase space factorisations for final-final, initial-final and initial- 
initial configurations have been derived in [32,79,82]. 

The presence of massive hard radiators introduces slight differences with respect to the 
massless case in the factorisation formulae for final-final and initial-final configurations. In 
the initial-initial configuration, the situation is unchanged as the radiators are massless. 
In the following, we shall present the factorisation formulae for the phase space in these 
three configurations while restricting ourselves to the kinematical situations where the two 
unresolved partons are massless. 

The final-final case has already been derived in [90] . We recall it here for completeness. 
We start from the 2 — )• (m + 2) parton phase space d$m+2 as given in eq.(2.3) involving 
two massive and m massless partons whose momenta are denoted by p3, . . . ,Pm+4- It can 
be written as 

d$m+2(P3, ■ • • ,Pm+4;Pl,P2) = d^rn{P3, ■ ■ ■ ,PI,PL, ■ ■ ■ , Pm+2; Pi , P2) 

X d^XijkliPi^Pj^Pk,Pi;Pl + Pl)- (2.12) 

In this case, pi and pi are the momenta of the massive radiators i, I between which the 
colour-connected and massless unresolved partons j, k with momenta andpk are emitted. 
The antenna phase space d^Xijki{PiiPjiPk^Pl\Pi + Pl) is closely related to the 4-particle 
massive phase space d^4{pi,pj,pk,pi;pi +Pl) as follows 

d'^4{Pi,Pj,Pk,Pi;Pi +Pl) = P2(.q^,mQ) x d^Xij^tiPi^Pj^Pk^Pl^Pi + Pl), (2.13) 
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where = {pi and P2{q'^,mQ) denotes the inclusive phase space measure for two 

particles of equal masses mg in d = 4 — 2e dimensions. 



For the subtraction terms built with massive initial-final four-parton antenna functions 
(with a massless initial state hard radiator and a massive final state one) we want our 
2 — )• (m + 2) phase space to factorise into a 2 — )• m reduced phase space (with massive and 
massles final state particles) times a 2 — t- 3 antenna phase space with at least one massive 
final state particle. To show how the factorisation is derived, we start again from the same 
2 — >■ (m + 2) particle phase space and consider the following situation: the massless initial 
state hard radiator has momentum pi , the massive final state radiator has momentum pi , 
the unresolved and colour-connected particles have momenta pj and pk- We insert 

1 = j dU5{q+pi-pj-pk-pi), (2.14) 
jW^l]{2T,f5{q + xpi-p^l) (2.15) 



and 



27r ./ X 



with 



and 



Q2 = -g2 ^2^^ = ^2 + ^2 ^ ^2 (2.16) 

(2.17) 



^ _ Q ^'^jkl _ Slj + -Sj/, + .Sj/ - .S;7, - Sji - Ski 

2pi • q sij + sik + sii 

in eq.(2.12), the expression of the 2 ^ (m -|- 2) phase space. By doing so, we find that it 
can be conveniently written as follows 

d$m+2(P3, • ■ ■ ,Pm+4;Pl,P2) = d$rn(Pl, • • • ,Pjkl, ■ ■ ■ ,Pm+2] Xpi,P2) 

X ^-^d$3(Pj,Pfc,P/;Pi,g) — • (2.18) 

ZTT X 

The four-parton intial-final phase space d^Xij^i involving one massless and one massive 
radiators in the intial and final state respectively, given for z = 1 by 



'^^x,,M = ^:^dMPj,Pk,Pi;Pi,q), (2.19) 



can be regarded as the new initial-final massive antenna phase space. It is worth noting 
that, in the case where j, k are massless and only I is massive, we have rrijki = mi = niL- 
Moreover, another important observation is that, when written in terms of Sab = '^Pa ■ Pb, 
the rescalling variable x defined in eq.(2.17) has the same form as in the massless case [80]. 
Since this same variable x is used in the definition of the phase space mappings for initial- 
final double unresolved configurations (see appendix B.2.2), the aforementioned observation 
suggests that the phase space mappings should be the same in the massless and massive 
cases provided that all variables are expressed in terms of invariants as in eq.(2.17). This 
is indeed the case. 



- 15 - 



For the initial-initial phase space factorisation, we can use the phase space factorisation 
derived in [79, 82] as in the massless case. The hard radiators situated in the initial state 
are obviously massless, and massless (and massive) final state partons have to be redefined 
using a Lorentz boost. In this case the phase space can be factorised as 

d^m+2{P3, ■ ■ ■ ,Pm+4;Pl,P2) = d^rniPs, ■ ■ ■ ,Pm+4; XlPl, X2P2) 

x5{xi - xi) S{x2 - £2) [dpj][dpk] dxi dx2 , (2.20) 
where the tilde over the final state momenta indicates that they are boosted. 

Using the factorised forms of the phase space in the three configurations as defined 
above in eqs. (2.12), (2.18) and (2.20) one can explicitly rewrite each of the four-particle 
subtraction terms present in da^^^^ in the forms, 

\Mm+2\' Jt^ d^m J d$X,,H X^Jkl, (2-21) 
\Mm+2\'j!T^ d<^mJ^^^^dMpj,Pk,pi;p,q) Xfj,i^, (2.22) 
\Mm+2\'' Jt^ d^m J [dpj][dpk]S{xi - Xi) S{X2 - X2) Xfi j^dXi dX2. (2.23) 

Using this, we then carry out the integration of the corresponding four-parton antenna 
functions. The integrated forms of the four-parton antennae are defined by 



1 f Q H~ ^ 

Kjkli^i) = ^2 J d^^^^^ljkl, (2-25) 



where 



C{e) = {AnY^, (2.27) 

and where the antennae and antenna phase spaces including possibly the masses of the 
radiators have been derived above in eqs. (2.12), (2.18) and (2.20). 

These integrations are performed analytically in d dimensions to make the infrared 
singularities explicit. 

2.4 Subtraction terms for colour-unconnected unresolved radiation duj^^j^Q 

We finally arrive at the last configuration considered in this section and relevant for our 
purpose in this paper, namely the configuration [d) involving the situation where the two 
unresolved partons are pair-wise colour-unconnected. For colour-ordered amplitudes of the 
type ..,i,j,k,...,n,o,p,...), where partons j and o arc unresolved, the (m+4)-parton 
matrix element factorises into the product of two uncorrclated single unresolved factors with 
hard partons /, K and N, P respectively multiplied by a reduced (m + 2)-parton matrix 
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element. The structure is the same as in the massless case [82]. The difference arises only 
from the use of three-parton antennae with massive radiators and massive phase spaces. 
The subtraction terms for the colour-unconnected configuration read, 



1 



perms 



'm+2 



^ XO., Xl,p \Mm+2{. ..,I,K,...,N,P,...)\^ 

3,0 

XJ^\PS, ...,PI,PK,--- ,PN,PP, ■ ■ ■ ,Pm+4) 
^^NNLO = --^NNLO d^m+2iP3, ■ ■ ■ ,Pm+4;Pl,P2)-^ 



perms 



E E ^i,jk <op \Mm+2i. ..,i,K,...,N,P,... 

i=l,2 j,o 



Sm+2 
2 



^4r^(P3, ...,PK,---,PN,PP,--- ,Pm+4) 



^^NNLO = --^NNLO ^ d*m+2 (P3 , ■ • • , Pm+4 ; Pi , P2 ) 

tJi — 



perms 



m+2 



J2 E ^M-fc Kop \Mm+2i. ..J,K,...,N,P,...)\' 

i,n=l,2 j,o 



'~jt\p^ 



,PK,---,PP,---,Pm+4) 



+ E J2^ljiKop\Mm+2{---,I,K,...,N,P,...)\'' 

k,n=l,2 j,o 

XJ^\P3, ...,PI,...,PP,... ,Pm+4) 

+ E J2 ^ik,J^nop\Mm+2i---J,K,...,N,P,...)f 

i,k=l,2 j,o 



xJ^Hp3, ...,Pn,PP,--- ,Pm+4) 



(2.28) 



(2.29) 



(2.30) 



where the summation over o is such that it only includes two antenna configurations with 
no common momenta. The nature of the radiator pairs {i, k) and {n,p) defines the antennae 
to be used. As it was already mentioned before, this term compensates spurious double 
unresolved features arising in subtraction terms for configurations (a) and for terms in (6) 
proportional to products of two three-parton antenna functions. 
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To obtain the integrated form of this counterterm we exploit the factorisation of the 
(m + 2)-parton phase space into an m-parton phase space multiphed by two independent 
phase space factors for each of the two antennae. The integrated form is thus simply the 
product of two integrated three-parton antennae. Finally, as discussed in [83] this subtrac- 
tion term has to be added back in integrated form at the m-parton level. 

As it was already the case at the NLO level [89], the subtraction terms derived above 
in all configurations entering at NNLO, are strictly valid for the subtraction of infrared 
singularities of colour-ordered matrix elements squared. Interferences between partial am- 
plitudes with different colour orderings appearing in the subleading colour pieces can also 
develop infrared singularities which need to be subtracted. As we will see in section 5.2, 
the subtraction of infrared singularities in those terms does require special care. However, 
to keep our equations as brief and clear as possible, we still write our subtraction terms in 
all configurations (f-f,i-f and i-i) and in all cases (a), (6) or (d) as above. 

3. Antenna functions 

In this section, we present the expressions of the four-parton antenna functions which 
are explicitly needed in the subtraction terms for the double real radiation corrections to 
hadronic heavy quark pair production due to the partonic processes obtained by crossing 
any pair of massless fermions in the process — t- QQqqq'q'- The two massless quark- 
antiquark pairs can (but do not necessarily have to) be of the same flavour. 

As explained in the previous sections, antenna functions are the key ingredients needed 
to build subtraction terms in the antenna formalism. In general, those functions are denoted 
with the character X. Each antenna is determined not only by its particle content (hard 
radiators and unresolved particles) but also by the pair of hard particles that it collapses 
to in its singular limits. Antennae that collapse onto a qiiark-antiquark pair are X = A for 
qgq and qggq, X = B for qq'q'q and X = C for qqqq. Antennae that collapse onto a (anti) 
quark and a gluon are X = D for qgg and qggg, and X = E for qc^q' and qq'q'g. Finally 
gluon-gluon antennae are X = F for ggg and gggg, X = G for gqq and gqqg, and X = H 
for qqq'q' ■ 

All antenna functions are derived from physical colour-ordered matrix elements squared. 
Each type of antenna is obtained from a different physical process: quark-antiquark an- 
tennae are related to processes of the form 7* — t- gg+(partons) [65], quark gluon antennae, 
to X ~^ 5+(partons) [66], and gluon-gluon antennae are derived from H — >-(partons) [67]. 
The tree-level antenna functions are obtained, both in the massless or massive case, by 
normalising the colour-ordered three- (or four)-parton tree-level squared matrix elements 
to the squared matrix element for the basic two-parton process, omitting couplings and 
colour factors: 




(3.1) 



(3.2) 
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S denotes the symmetry factor associated with the antenna, which accounts both for po- 
tential identical particle symmetries and for the presence of more than one antenna in the 
basic two-parton process. 

The three and four-parton initial-final and initial-initial antenna functions are, in prin- 
ciple, defined by crossing one or two massless partons from the final to the initial state 
in the corresponding final- final antennae. The three-parton initial- final and initial-initial 
antennae are denoted respectively by X^jj. &iid X^j^j while the corresponding four-parton 
antennae will be denoted by X^^j^i and X?^ ^ . 

While at NLO only three-parton tree-level antennae are used, at NNLO also four- 
parton tree-level and three-parton one-loop antenna functions need to be considered. The 
latter, however, are not needed for the double real radiation contributions, and they will 
not be discussed here. All the three-parton antenna functions used in the subtraction terms 
that will be presented in section 5 are of types A and E in different configurations with 
massive and/or massless partons. These antennae have been computed, integrated and 
extensively studied in [32,79,88,89], and we shall therefore not discuss them further here. 
They are listed in appendix A together with their soft and collinear limits. 

We shall here focus only on the four-parton antennae which are genuine NNLO objects 
capturing double unresolved radiation features of real matrix clement squared. In the 
context of this paper, the double unresolved limits which need to be subtracted with 
four-parton antennae are double soft and triple collinear limits involving only fermions in 
initial and final states. All these limits can be subtracted with B-type and C-type antenna 
functions. The former are needed in their final-final, initial-final and initial-initial forms, 
while only initial-initial massless forms of the latter are needed. The B and C-type antennae 
only exist at the four-parton level. As the A-typc antenna functions, and antennae 
collapse onto a hard quark-antiquark pair in their double unresolved limits, but they have 
a secondary qq pair in the final state instead of gluons. More precisely, in the final-final 
case, B2 and antennae are derived from 7* — > qqq'q' [32]. While is obtained from 
a colour-stripped amplitude squared, C4 is a pure interference term. Explicitly, they are 
given by 

b2(1„V,3,,,2,-)|A^o/ = |7W0(lg,V,3,.,2,-)|', (3.3) 
C0(l„3„4,-,2,-) |A^o/ = -Re (A^O(lg, 3^., 2,-)A^0(l„ 2,-,, 3,., 4,-)) , (3.4) 

where the amplitude A4^{iq,lqi , kg/ , jq) contains those Feynman diagrams contributing to 
the colour stripped amplitude for 7* — t- iqjqkqdqi in which the pair iqjq is the one that 
couples to the virtual photon. Explicit expressions for the final-final massless antennae B^ 
and C4 are given in [32]. 

As mentioned before, the initial-final and intial-initial antennae are obtained by cross- 
ing partons in the final-final functions. Their infrared limits are in principle known as they 
can also be inferred from the final-final case, although they have not been documented so 

^It is woth noting that, in the initial- final case this crossing may not be unambiguous and may result 
in the necessity of a further splitting of the antenna as noted in [79,89]. However, those antennae are not 
needed in the present context of this paper. 
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far. Furthermore, concerning the four-parton antennae involving massive partons, only the 
final-final massive is known [90]. 

In addition to the massive flavour conserving A-type (involving three-parton final 
states) or B-type antennae (involving four-parton final states) which can be derived from 
the decay of a virtual photon, the construction of our subtraction terms will also require the 
so-called massive flavour violating antennae. These involve radiators of different flavours, 
one of them being massive and the other massless. For the three-parton case, flavour vi- 
olating antennae have been derived and integrated in [89]. The required massive flavour 
violating antenna with four final state partons is new and will be given below. 

In the remainder of this section we give the forms of the four-parton antennae used in 
our calculation. Their infrared limits will be presented in section 4. The conventions that 
we will follow for the labelling of the partons in our antenna functions are the following: 
Massless partons will be indexed with q while massive ones with Q and their mass with 
ruQ. The first and the last particles in the argument of a given antenna are the hard 
radiators, while the partons placed in between the radiators are the unresolved particles. 
For conciseness, the 0{e) in the expression of the unintegrated antennae will be omitted. 

B-type antennae 

The subtraction of the double unresolved singularities in the processes that we are presently 
considering requires B-type antenna functions in their massive final-final, fiavour violating 
initial-final, and massless initial-initial forms. They are given below. 

Final-final 

The only final-final antenna function which is needed in our subtraction terms is the 
one involving a pair of massive radiators (QQ) and a secondary massless quark-antiquark 
pair (qq). It is denoted by i?4(lQ, 4^, 3^, 2q), and it has been computed and integrated 
in [90] . This antenna function is used to subtract the infrared limits in which a soft quark- 
antiquark pair is emitted between a massive quark and a massive antiquark in the colour 
chain. Its explicit expression reads, 

-B4 (1q, 4g, 3g, 2n) = 7 ^ j (si2Si3 + si2Si4 + S13S23 + S14S24) 

H K— (S12S23 + S12S24 + S13S23 + S14S24) 

S34S234 

H — 2 — 2~ (2S12S13S14 + S13S14S24 + S13S14S23 - sl^S24 " Sl4'S23) 
*34*134 

+ ^ 2^ (2S12S23S24 + S13S23S24 + S14S23S24 - ■Sl3si4 - S14S23) 

H (2si2 + S12S23 + S12S24 + S12S13 + S12S14) 

S34S134S234 

+ ^ ( - 513514 - •S14SI3 - 513524 - 514S23 

5§4S1345234 
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+S13S14S23 + S13S14S24 + S13S23S24 + S14S23S24 



S134S234 



-2si2Si3S24 - 2S12S14S23) + 

2 f 8S13S14 8523^24 4 ^ 

+"^<3l „2 „2 + „2 „2 „2 „2 

\*34*134 *34*234 *134 *234 

2 2 

2—{si2 + S23 + S24) —{^12 + Sl3 + Su) 

S34S134 ■534S234 

2 

H (4si2 - Si3 - Su - S23 - S24) 

S34S134S234 

8 

(S14S23 + S13S24) 

S34S134S234 



-<J 2- + 2- + ^(^)' (3-5) 

Vs34Sf34 5345^34/ J 

with the normaUsation given by the tree-level two-parton matrix element (with couplings 
and colour factors omitted) 

\Ml{j* QQ)\' = 4 [(1 - 6)£;2^ + 2ml] ■ (3-6) 

Initial-final 

To account for the singularities that arise when a massless quark-antiquark pair becomes 
unresolved between a massive final state fermion and a massless incoming one, we need 

a massive initial-final flavour-violating B-typc antenna. The hard radiators are a massive 
final state quark Q, and a massless initial state q, and the expression of this antenna is 

-64(1(5, 4q,3q,2q) = --, ^ (si2Sl3 + Sl2Sl4 + 513^23 + S14S24) 

[Q^ + mlj I ^34^134 

H ^^o— (S12S23 + S12S24 + S13S23 + S14S24) 

S34S234 

- 22 (2'Sl2Sl3Sl4 + S13S14S24 + S13S14S23 - S13S24 - S14S23) 
*34*134 

H 2 2 (-2S12S23S24 + S13-523S24 + S14S23S24 - S13SI4 - Susls) 

•^34*234 

H (25^2 + S12S23 + S12S24 - S12S13 - S12S14) 

S34S134S234 

H 2 ( ~ ■513S24 - ■S14S23 + S13S24 + S14S23 

S34SI34S234 

-S13S14S23 - S13S14S24 + S13S23S24 + Sl4S23'S24 

\ 2si2 

-2si2Sl3S24 - 2S12S14S23) 

S134S234 

+2ml ( -i- + + '-11^) ] + 0{e), (3.7) 

\'Sl34 S134S234 S34S134 S34Sf34 yj 
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where = -(pi - P2 + P3 + Pa)"^, •S134 = si3 + su + S34, and S234 = -S23 - S24 + S34- In 
this case the normahsation is 

\M^2{l*Q ^ = 4(1 - e)(Q2 + ml). (3.8) 



Initial-initial 

Two initial-initial forms of B-type antennae are needed in our subtraction terms. The 
first of them is -64(19, 4^/, 3g/, 2g) and it used to capture the limits where a final state quark- 
antiquark pair becomes unresolved between the incoming quark and antiquark, which act 
as the hard radiators. As will be shown in section 4, this antenna contains the 3^/, 4^/ double 
soft limit as well as the triple collinear limits lq||4q/||3q/ and 4q/||39/||2g. The expression of 
this first initial-initial antenna function is, 

B2{iq,4q/,3q',2q) = -"^j (*12Sl3 + ■S12S14 + S13S23 + S14S24) 

V L S34S134 

1 / 

2— (S12S23 + S12S24 + S13S23 + S14S24) 

*34'S234 

(2si2Sl3Sl4 - S13S14S24 - S13S14S23 + Sl3'S24 + S14S23) 



SoaS 



34 •'134 



H 2 ^2 (2S12S23S24 - S13S23S24 - S14S23S24 + Sissli + S14SI3) 



34 •'234 



H (2Si2 - S12S23 - S12S24 " S12S13 - S12S14) 

S34S134S234 

H 2 (S13S24 + «14S23 + Sl3«24 + S14S23 

•S34'Sl34S234 

-S13S14S23 - S13S14S24 - Sl3S23'S24 - S14S23S24 

-2si2Si3S24 - 2S12S14S23) H ^^^^ \ + C(e), (3.9) 

S134S234 J 

with = -{pi + P2 - P3 - Pa)"^ , •S134 = -Si3 - + S34, and S234 = -S23 - ^24 + S34- This 
antenna function was computed and integrated in [81]. It is normalised to 

\Ml{qq^l*)\'' = 4{l-e)Q\ (3.10) 

The second B-type initial-initial antenna function that we use in our subtraction terms is 
B2{lq, 4:qi,3qi, 2q) and it only contains the lq||4q/||3g/ triple collinear limit. The unintegrated 
form of this antenna is 

B4{lq,4:qi,3q',2q) = (si2Sl3 - Sl2Sl4 + S13S23 + S14S24) 

V LS34S134 

1 / 

-I 2~ (-•S12S23 + S12S24 + S13S23 - S14S24) 

^34*234 

H 2 ^2 (2S12S13S14 - S13S14S24 + S13S14S23 - S13S24 + Sl4'523) 

•*34^*134 
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„2 2 
* 34 -^234 



(2si2S23S24 - S13S23S24 + S14S23S24 - S13SI4 + Suslg) 



(2si2 + S12S23 - S12S24 - S12S13 + S12S14) 

S34S134S234 

H 2 ( ~ *13S24 + SUsis - S?3S24 + S?4S23 

•S34'Sl34'S234 

+S13S14S23 - S13S14S24 - S13S23S24 + S14S23S24 

+2S12S13S24 + 2S12S14S23) - \ + 0{e) (3.11) 

S134S234 J 

with Q2 = -{pi - P2 + P3 - Pif, S134 = .Sl3 - Sl4 - ^34, and .S234 = -S23 + ^24 - S34- 

Its integrated form has been given in [81], and the normahsation is the same one given in 
eq.(3.10), but with the corresponding to this case. 

C-type antennae 

To subtract triple collinear limits in amplitudes with two qq pairs of the same flavour, two 
different initial-initial C-type antenna functions involving two massless final and two initial 

state partons are needed. 

The first of these initial-initial C-type antennae is 

/^O/? o o 3 ^ / S12S13S14 , 1 / , X 

Cl(lq,\,2q,Aq) = — ^ <^ V- (Sl2Sl3S24 + S13S14S24) 

y S23S34S123S134 2S23S34S134S234 
+7, (S12S13S24 + S13S14S24J 

2S23'S34Sl23S234 

Sl3'S24 _^ S12S13 S13S14 



S23S34S234 S23S123S134 S34S123S134 

+ 7. (■S12S14 - S12S34 - S12 + S13S24) 

2S23S123S234 

+ 7. (S12S14 - S14S23 - + S13S24) 

2S34S134S234 

(-S12S14 + SI2S34 + S12 + SI3S24) 

2S23S134S234 

-7. (-S12S14 + S14S23 + -I- S13S24) 

2S34S123S234 

si3 , 1 / ^ 1 ^ ^ 

\ 2~ (*12S24 - Sl4S24j o" (*12S24 - ■Sl4S24j 

Sl23'Sl34 ■S23S234 -5345234 

+ 7r-^ (si2 + S14) + ^ (si2 + S14) j + 0(e), (3.12) 

^S123S234 ^Sl34S234 J 



which is used to subtract the triple collinear limit 3g||2g||4^. The second one is given by 

r^Of^ 6 o /I ^ ^ / si2Si3'Si4 1 , s 

W(lg)^g)25,45j = (S12S13S24 - S13S14S24) 

Q y S23S34S123S134 2S23S34S134S234 



1 

(-•Sl2Sl3'S24 + Sl3Si4S24j 



2S23S34S123S234 

Sl3si4 S12S13 S13S14 



S23S34'S234 S23Sl23Sl34 S34Sl23Sl34 
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+7, ^ (-S12S14 - S12S34 - S12 + S13S24) 

2S23S123S234 

+ 7, (-S12S14 - S14S23 - + S13S24) 

2S34S134S234 

+ -Z (si2SU + SI2S34 + S12 + SI3S24) 

2S23'S134S234 

+ 7, (S12S14 + S14S23 + •S14 + S13S24) 

2S34S123S234 

■S13 1 / ^ 1 / \ 

\ 2~ (*12S24 + Sl4S24j H (*12S24 + Sl4S24j 

■Sl23'Sl34 ■S23S234 •S34'S234 

+ (-S12 + S14) + (si2 - su) ] + 0{e), (3.13) 

^Sl23S234 ^^1345234 J 

and it is used to subtract the triple coUincar limit 3g||2q||4q. 

The normalisation for the C-typc initial-initial antennae is again given in eq.(3.10) with 
given by the square of the sum of incoming parton momenta in each case. 



4. Infrared limits 

The factorisation properties of QCD tree-level squared amplitudes have been extensively 
studied in [93 -96]. While at NLO only single soft and collinear singularities may arise, at 
NNLO there are several double unresolved regions of phase space in which the amplitudes 
can develop singularities when integrated over. 

4.1 Single unresolved factors 

In their infrared limits, colour-ordered matrix elements factorise into universal unresolved 
factors and reduced matrix elements. Depending on whether the partons involved in each 
limit are massive or massless, and depending also on whether one or two partons become 
unresolved, the universal factors will be different. We shall start by presenting the single 
unresolved massless and massive factors. 

Retaining only those infrared limits which lead to poles in e when integrated over the 
phase space, the only single unresolved limits that colour-ordered amplitudes involving 
only fermions have are single collinear massless limits arising when a quark-antiquark pair 
becomes collinear and clusters to form a gluon. If both collinear particles are in the final 
state, the parent gluon will also be in the final state, and it will be in the initial state 
in the cases where one of the collinear partons is in the initial state. Thus, the reduced 
colour-ordered matrix element associated with this limit involves a gluon in the final or 
in the initial state. Therefore, although only quarks are involved in the matrix elements 
for the processes that we are considering in this paper, unresolved gluon limits have to be 
considered as well. While the collinear limits of this gluon only give poles in e if the other 
collinear particle is massless (a massless fermion in this case) , in soft limits mass terms will 
appear in the soft eikonal factors if the gluon is emitted from a massive hard radiator. 
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4.1.1 Single unresolved massless factors 

When two massless partons become collinear, a colour-ordered sub-amplitude factorises into 
a reduced matrix element times a specific Altarelli-Parisi splitting function corresponding 
to the particular parton-parton splitting. These functions depend on z, the momentum 
fraction carried by the unresolved parton. Depending whether the unresolved parton is 
collinear to an initial or to a final state parton, the definition of z is different. For two final 
state particles i and j of momenta and becoming collinear, we have, in the limit, 

Vi zpij, Pj (1 - z)pij, Sik zsijk, Sjk (1 - z)sijk , (4.1) 

whereas for a final state particle j of momentum pj becoming collinear with an initial state 
parton i of momentum pi we have 

Pj ^ zpi, pij^{l-z)pi, Sik^Y^, ^jk^YZ^' ^^-'^^ 

The splitting functions denoted by P,ij^(^ij^ [z) corresponding to the collinear limit of two 
final state partons i and j, given in the conventional dimensional regularisation scheme [97] 
with all particles treated in d = 4 — 2e dimensions, are [91]: 

(4.3) 

P.-.a(.) = (4.4) 



Pgg^Giz) = 2 



z 1 — z . 
+ + z{l - z) 



(4.5) 



When the collinearity arises between an initial i and a final state parton j, the splitting 
functions denoted by Pij^(^ij^{z) are given by [91]: 

, , 1 + z'^ -e(l- zf 1 1 „ , , , 
W^)= = 13^131 Wl-^) (4.6) 

1 + (1 - z)^ - ez^ 1 

P,9^Qi^) = \^ = Y^/,9^Qi^) (4.7) 

z"^ -\- (I — zS^ — f 1 — e 
P,,^g{z) = ^\_"J = —P,,-,g{z) (4.8) 

2(1-2 + ^2^)2 ^ 

Pgg^Giz) = \^^_^^2 = Y^/99-^g{z). (4.9) 

The additional factors (1 — e) and 1/(1 — e) account for the different number of polarizations 
of quark and gluons in the cases in which the particle entering the hard processes changes 
its type. 

In all the splitting functions defined above, the label q can stand for a massless quark or 
an antiquark since charge conjugation implies that Pqg^q = Pqg^Q and Pqgi^Q = Pgg^Q ■ 
The labels Q and G denote the parent parton of the two collinear partons, which is massless. 
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When a gluon j emitted between two masslcss hard radiators i and k becomes soft, the 
eikonal factor that factorises off the colour-ordered squared matrix element is 



Sijk = (4.10) 

This limit is obtained by letting pj Xpj with A ^ or, alternatively, by keeping only 
terms with any product of the inverse of the two invariants Sij and Sjk in the colour-ordered 
matrix element squared. 



4.1.2 Single unresolved massive factors 

In the present context, the only single unresolved massive factor that needs to be considered 
is the massive soft factor, which is a generalisation of the massless eikonal factor defined 
above. The generalised soft eikonal factor Sijk{mi,mk) for a massless gluon j emitted 
between two massive partons i and k depends on the invariants sim = "^Pi ■ Pm built with 
the partons i, j and k but also on the masses rui and of i and k. It is given by [85,88] 

Sijk{mi,mk) = 2^ (4.11) 

Sij Sjk S-j Sjf^ 

This limit is again obtained by letting pj — > Xpj with A — t- or alternatively by keeping in 
the matrix clement squared only terms with the inverse of any product of the invariants 
Sij and Sjk- The above expression for the massive soft factor agrees with eq.(4.10) if the 
masses are set to zero. 



4.2 Double unresolved factors 

In general, when two particles are unresolved in a tree-level process, a variety of different 
configurations can arise: 



(1) two soft particles, 

(2) two pairs of collinear particles, 

(3) three collinear particles, 

(4) one soft and two collinear. 



As we saw in section 2, the resulting double unresolved configurations of the colour- 
ordered matrix element squared need to be separated into three categories depending on 
the colour connections of the unresolved partons and the hard radiators associated with 
them. In the following, we shall describe only the double unresolved factors encountered 
in the context of this paper. 
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4.2.1 Colour-unconnected pairs of unresolved particles 

For unresolved particles that are disjoint in the colour chain, which arise in the item (d) 
as presented in section 2, the colour-ordered matrix element can be factorised into the 
product of two disjoint single unresolved factors multiplied by a reduced matrix element 
with two partons less than the original colour-ordered matrix element squared. 

For example, in the case of two pairs of colour-connected particles (a, 6) and (c, d) 
which are all located in the final state, the colour-ordered matrix elements given by 
(. . . , a, 6, . . . , c, d, . . .)p undergo the following factorisation in the double collinear limit: 

(. . . , a, 6, . . . , C, d, . . .)P ^ Pab^p{zi, Sab) Pcd^Q{z2, Scd)\Ml_2{. . . , P, . . . , Q, . . 

(4.12) 

where partons a and b form P, while c and d cluster to form Q, so that P and Q are 
themselves colour-unconnected. The collinear splitting functions appearing in this equation 
are the (colourless) Altarelli-Parisi splitting functions given by eqs. (4.3-4.9). 

In the processes that we are considering in this paper, the only double unresolved 
colour-unconnected limits that can occur are double quark-antiquark collinear limits, with 
each single collinear pair given by a final state (anti) quark and an initial state (anti) quark. 
The unresolved factor associated to this limit is a product of two splitting functions of the 
type given in eqs. (4.3-4.9). Note that none of the four-parton antennae required for our 
subtraction term captures these double collinear singularities. These are accounted for by 
subtraction terms involving only the product of two three-parton antennae. 



4.2.2 Colour-connected pairs of unresolved particles 

The colour-connected double unresolved limits that appear in the partonic processes con- 
sidered here are 

A) Triple collinear limits of three massless fermions where one massless fermion plays 
the role of the hard radiator and is in the initial state. 

B) Double soft limits of a massless final state qq pair emitted between massive or massless 
radiators. 

As it was mentioned above, these colour-connected double unresolved limits are captured 
by four-parton antenna functions and arise in subtraction terms dcr^^^^ as presented in 
section 2. In the following, we shall give the massive and massles double unresolved factors 
associated with these two types of limits. 

A) Massless triple collinear factors 

In those regions of phase space where three colour-connected massless partons (a, 6, c) be- 
come collinear, a generic colour-ordered amplitude squared denoted by |A^° (. . . , a, 6, c, . . .)p 
factorises as: 

(. . . , a, 6, c, . . .)|2 ^ Pabc^p\M^n-2i- ■■,P, ■■■){''■ (4.13) 
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where the three colour-connected final state particles (a, b, c) cluster to form a single parent 
particle P. The triple collinear splitting function for partons a, b and c clustering to form 
the parent parton P is generically denoted by, 

Pabc^pi'U}, X, y, Sab, Sac, Sbc, Sabc), (4-14) 

where w, x and y are the momentum fractions of the clustered partons, 

Pa = wpp, pb = xpp, Pc = ypp, with w + x + y = l. (4.15) 

In addition to its dependence on the momentum fractions carried by the clustering partons, 
the splitting function also depends on the invariant masses of parton-parton pairs and the 
invariant mass of the whole cluster. The explicit forms of the triple collinear splitting 
functions Pabc^p are obtained by retaining terms in the colour-ordered matrix element 
squared that possess two of the 'small' denominators Sab, Sac, Sbc and Sabc- 

The triple collinear limits that have to be considered here are those involving three 
massless fermions, one of which is in the initial state. The splitting functions for these 
types of triple collinear limits can be obtained from the analoguous limit where the three 
collinear particles are in the final state [93,96]. 

We therefore start by considering the case where a quark-antiquark pair q'^ and a 
quark q of different flavour cluster together to form a parent quark q of the same flavour as 
q. In such a configuration a squared sub-amplitude in which the three collinear fermions 
are colour-connected factorises as: 

(4.16) 

with pi = Pi + pj +pk- The corresponding splitting function is 



^aa'a' V / 



qq'q 

2{xSqq>g> - {l-w))Sqq>)'^ ^ 1 / l + x'^ + {x + wf _ _ 

(4.17) 

where w, x, and y are the momentum fractions of q, q' and q' respectively, and w+x+y = 1. 

As we will see in section 5, the squared amplitudes for processes with identical quark 
flavours contain interferences of partial amplitudes with different flavour assignments for 
the massless antiquarks. Although they are finite in all other single and double unresolved 
limits, these interference terms can still be singular in triple collinear configurations. In 
these cases, the factorisation reads 

iq\\jqi\kq 



Re (M'l{...,ig,jg>,kg>,m,...)Mni---,iq',jq',kq,n,...)^^ 

-^PqPq^^qiw, X, y, Sqq> , Sg/g> , Sgq>g/)Re (^A^°_2("-, Iq, "^, ...)A^°_2(..., Iq, U, ... 



(4.18) 
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The corresponding splitting function, given by 

1 



^gg'gA V?' / "^^qq' Sqi \{l - y){l - w) \l-W 

+ e ^- + 1 + X + + e(l - w, 



Sq'q'Sqq'q' \ l-y l-W \ l-y l-W 

+ {Sqq> O Sqiqi,y^ w), (4.19) 

is manifestly symetric under the interchange q q' and y w. 

The triple coUinear splitting functions in eqs. (4.17) and (4.19) correspond to configura- 
tions in which all three coUinear particles are outgoing. However, having only two massless 
final state particles, the double real radiation matrix elements for top-antitop production 
can only have triple coUinear limits involving one of the incoming partons and both mass- 
less outgoing particles. The corresponding initial-final splitting functions are related to 
their final-final counterparts by [96] 

Pqq' q' *^q{,Z\, ^2) ^3) ^qq' t ^q'q'i ^qq'q') — 

Pqq'q'^q{^/z3,—Zl/z3,—Z2/zs,—Sqqi,Sqiqi,Sq/q/ - Sqqi - Sqqi) 
Pqqq'^q'izi,Z2, Z3, Sgg' , Sg'g' , Sgg'g' ) = 

Pqq'q'^qi-Zl/Zs, I/Z3, -Z2/Z3, -Sqqi, -Sqfqi,Sqqi - Sqiqf - Sqqi), 

(4.20) 

where zi and Z2 are the momentum fractions of the final state particles participating in 
the limit, and z^ = \ — zi — Z2- eq.(4.20) is valid for both 'non-identical' and 'identical' 
triple coUinear splitting functions. 

B) Massive and massless double soft factors 

When a massless quark-antiquark pair becomes soft between two hard radiators, the am- 
plitude squared factorises into a double soft factor and a reduced matrix element squared 
with the quark-antiquark pair removed from it. As it was the case for the single soft gluon 
eikonal factor, the soft factor for a soft qq pair also depends on the masses of the hard 
radiators if these are massive. Again, the massless factor can be obtained from the massive 
one by setting the masses to zero. 

When the quark-antiquark pair (c, <£) becomes soft between the hard radiators (a, 6) of 
masses nta and m^, the massive double soft factor is given by. 



o / X 2{SahS(,d - SacSbd " SbcSad) , 2 



SacSad ^ SbcSbd 



{Sac + SadY {sbc + SbdY^ 



2m| 

Scd{Sac + SadY Scd{sbc + SbdY 

This factor is obtained by setting pc — )■ Xpc, Pd ~^ ^Pd with A — t- in the matrix elements 
squared and was also derived in [90]. This limit is only present in the processes qq — )• QQq'q' 
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and qq — t- QQqq. It occurs in those diagrams in which the final state quark-antiquark pair 
sphts from a gluon propagator that becomes soft. As it wiU be seen in section 5, we 
encounter all possible combinations of hard radiators (two massive, one massive and one 
massless, and two massless), and therefore eikonal factors with two masses, with one mass, 
and with no masses are needed. They can all be obtained from eq.(4.21). 

Furthermore, a derivation of this double soft factor using the soft current formalism 
will be given in section 5. 

4.3 Single and double unresolved limits of four-parton antennae 

In this section we present the single and double unresolved limits of all four-parton antenna 
functions defined in section 3. Those are the ones required to construct the subtraction 
terms for the double real contributions to heavy quark pair production in hadronic collisions 
due to processes involving only fcrmions. In their infrared limits, these antennae yield the 
universal single and double unresolved factors defined above. Additionally, some of these 
four-parton antennae yield angular correlation terms in single collinear limits. These an- 
gular terms arise when a gluon splits into a quark-antiquark pair or into two gluons. The 
treatment of these angular dependent terms will be discussed in section 4.4. 



The massive final-final B-type antenna given in eq.(3.5) is only singular when the mass- 
less qq is soft or when this quark-antiquark pair is collinear. The behaviour of this antenna 
in each of these limits is given by 



where (ang.) indicates the presence of angular correlation terms. 
Initial-final 

The flavour violating B-type antenna function given in eq.(3.7) contains a triple collinear 
singularity, in addition to the double soft and single collinear limits. It reads. 



Final-final 




(4.22) 



(4.23) 




(4.24) 




jnon— ident. 

qq'q'->^q 




(4.25) 




(4.26) 



where (ang.) means that angular terms are present here. 



Initial-initial 
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The initial-initial B-typc antenna with the massless hard radiators in the initial state 
has a double soft singularity, as well as two triple collinear and a single collinear limit: 

B^iU, ^',Sg',2g) '"^^ '5i342(0, 0) (4.27) 

V,3,s2,) P°^^r4^'^*-(Z1, ^2, ^3, 514,534,^134) (4.28) 

Bliiq, 4,-,, 3,,, 2,) P^^^i^t'^^^^ ^3, 524, 534, 5234) (4.29) 

Bl{iq, 4,-,, 3,,, 2,) —P,g^G{z)Al{i^, (34)g, 2,) + ang.. (4.30) 

534 

Again we see that the splitting of a gluon into a collinear quark-antiquark pair yields 
angular terms. 

The other initial-initial B-type antenna function does not have a double soft singularity, 
since one of the unresolved particles is crossed to the initial state. It only contains a triple 
collinear as well as a single collinear limit: 

i?0(lg, V, 3,,, 2,) ^2, ^3, 514, 534, 5134) (4.31) 

b2(1„ 4,,, 3,,, 2,) —Pgg^Q{z) + ang. (4.32) 

534 

As mentioned in section 3, the C-type antennae given in eqs. (3.12, 3.13) are used 
to subtract triple collinear limits in processes with two pairs of identical flavour massless 
quarks. These are actually the only singularities that these antenna functions contain. 
For the present calculation, we only need to consider initial-initial C-type antennae, whose 
limiting behaviour is, 

C4(lg, 35,2^,4^) > -Pqq?qr^_q{zi,Z2,Z3, 824,834,8234), (4.33) 

and 

C2{iq, 3q, 2g, Ag) lP^tl'tq'{zi,Z2, Z3, 523, 524, 5234)- (4.34) 

4.4 Angular terms 

As seen above, in single collinear limits where a gluon splits into a quark-antiquark pair 
(or into two gluons), antenna functions do not just yield a product of a spin averaged 
splitting function (4.3-4.9) and a spin averaged reduced matrix element: angular correla- 
tion terms are also present. This is a general feature of QCD amplitudes squared and it 
holds both for four-parton antenna functions and for the real radiation matrix elements. 
In these types of single collinear limits, arising from gluon splittings, neither the real radi- 
ation matrix elements squared nor the four-parton antenna functions yield the unpolarised 
splitting functions multiplied by spin-averaged reduced matrix elements. Instead, ampli- 
tudes squared (or, similarly, four-parton antennae) factorise into spin-dependent tensorial 
splitting functions contracted with a tensorial reduced matrix element [21,95]: 

\M'i{...,a,;;b„...)f ^ ^P,^^^Gi^,k^)\MU{--.c,, (4.35) 



-31 - 



In a slight abuse of notation, |A^J^_]^(..., c^, denotes the reduced matrix element 

squared without the average over the polarisations of gluon c. The spin dependent splitting 
function is given by 

ki-) = -g'" + 4^(1 - (4-36) 
and the coUinear limit is parametrised as 

pf = (l-.)p^-*i-^^. (4.37) 

n is a reference light-like four-vector needed to uniquely determine the transverse com- 
ponent k±, which satisfies k± ■ p=0 and k± ■ n=0. The coUinear limit is approached as 
k± 0, and, as it can be seen from eqs. (4.35) and (4.36), the angular correlation terms 
are proportional to 

^'^\MU{...,c„...)\l^. (4.38) 

Sab "^_|_ 

Since in the antenna formalism single unresolved limits of the amplitudes squared are 
subtracted with spin-averaged three-parton antenna functions multiplied with reduced ma- 
trix elements which do not contain the azimuthal terms produced by the spin correlations 
of the matrix element, it is clear that the subtraction terms do not have the same pointwise 
singular behaviour as the amplitudes squared. 

In the same line of thought, a similar subtlety is found in the subtraction terms for two 
colour-connected unresolved partons. As explained in section 2 these subtraction terms 
(dcJ^^^Q) are constructed as a difference of four-parton antenna functions and products 
of two three-parton antennae multiplied with reduced matrix elements squared. The for- 
mer are intended to subtract all singularities in the double unresolved region, while the 
latter ensure that the whole subtraction term is free of singularities in all single unresolved 
regions. If the four-parton antenna in eqs. (2.9-2.11) has single unresolved limits where an- 
gular terms are present, the product of the two three-parton antenna functions in the same 
equation will no longer subtract its single coUinear behaviour locally. In both cases (real ra- 
diation matrix elements and unintegrated subtraction terms), however, angular correlation 
terms are averaged out once the integration over the azimuthal angle between the momenta 
Pi {i = a,b) and p (the coUinear direction), has been carried out for all coUinear limits. 
Moreover, it can be shown [68,82,92] that in gluon splittings the functional dependence of 
the angular terms in eq.(4.39) on the azimuthal angle (j) is 

1 k'^k'^ 

' '-|A^o_i(...,Cg,...)|2 ^cos(2<^ + a). (4.39) 



7.2 K-'n-iV"-! ••Vl/ii/ 
■'ao "'± 

This suggests that spin correlations can be averaged out by combining two phase space 
points whose azimuthal angles are 

(I) and 0+ ^- (4.40) 
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It is only after this combination of phase space points that the subtraction terms are a 
good approximation of the double real radiation corrections in all its single and double 
unresolved limits. This procedure, which has proven to be extremely successful for the 
double real radiation corrections to gluon scattering [82], is the one that we shall follow 
in section 6 to show that our subtraction terms correctly converge to the real radiation 
matrix elements squared in all their unresolved limits. 

It is also important to note that angular terms in the four-parton antenna functions 
vanish once they are integrated over the corresponding antenna phase space. This fact 
implies that the integrated antenna functions do not have any azimuthal correlations. This 
is crucial for the cancellation of infrared poles when the subtraction terms involving these 
antennae have to be added back in integrated forms at the m-parton level. 

5. Top quark pair production at the LHC 

In this section we shall present the colour-ordered double real emission contributions to 
tt production at the LHC due to processes involving only fermions. Together with these, 
we shall give their corresponding antenna subtraction terms, which capture all single and 
double unresolved limits of the respective real radiation matrix elements squared. 

5.1 Structure of double real radiation contributions 

In a previous paper [89] we have presented the colour-ordered form of the real radiation 
corrections as well as the corresponding NLO subtraction terms for all the partonic pro- 
cesses contributing to tt+jct production at hadron colliders. Both the NLO real radiation 
corrections and the subtraction terms for tt+jet are essential ingredients for the study of 
the double real radiation corrections to heavy quark pair production at NNLO: on one 
hand, the matrix elements involved and their colour decomposition are the same, and, on 
the other hand, the aforementioned NLO subtraction terms already capture all the sin- 
gle unresolved infrared behaviour of the matrix elements, which is needed at NNLO in 
the d&^'^j^Q piece of the subtraction term. Indeed, at partonic level, the forms derived 
in [89] for the real radiation matrix elements together with their related subtraction terms 
for single unresolved limits can be taken over to the present paper without any changes. 
While the real radiation matrix elements are exactly the same, the NLO subtraction terms 
for tt-|-jet production have to be supplemented with additional genuine NNLO pieces in 
order to obtain subtraction terms that account for the double unresolved features of these 
six-particle amplitudes squared. 

We shall start by discussing the contributions of purely fermionic processes to the double 
real radiation corrections to tt production at the LHC, and by recalling from [89] some of 
the properties of the corresponding real radiation matrix elements. 

5.1.1 Generalities 

To facilitate the reading of our expressions, we shall closely follow the notation in [89] 
for matrix elements and subtractions terms. The main points of our conventions are the 
following: The matrix elements denoted with Ai represent colour-ordered sub-amplitudes 
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in which the coupling constants and colour factors are omitted. Furthermore, to explicitly 
visualise the colour connection between particles in these colour-ordered amplitudes, a 
double semicolon is used in the labeling of the partons present in a given matrix element. 
This double semicolon is used for separating chains of colour-connected partons. Partons 
within a pair of double semicolons belong to a same colour chain, and adjacent partons 
within a colour chain are colour-connected. An antiquark (or an initial state quark) at the 
end of a colour chain and a like flavour quark (or initial state antiquark) at the beginning of 
a different colour chain are also colour-connected since the two chains merge in the collinear 
limit, where the qq is clusters into a gluon. Since in the antenna framework a parton can 
only be unresolved with respect to its colour-connected neighbours, this notation helps to 
identify the unresolved limits present in a given colour-ordered amplitude, and therefore 
helps to construct the corresponding subtraction terms. 

Notationwisc, we also denote gluons which are photon-like and only couple to quark 
lines, with the index 7 instead of g. In sub-amplitudes where all gluons are photon-like 
no semicolons are used, since the concept of colour connection in not meaningful in those 
configurations. A hat over the label of a certain parton indicates that it is an initial state 
particle (for example, 1,^ is an initial state quark with momentum pi). 

Concerning the notation in the subtraction terms, the conventions for the reduced ma- 
trix elements arc the same as those for the real radiation matrix elements discussed above. 
The remapped final-state momenta are denoted with tildes and the remapped momenta of 
initial state hard radiators are denoted by a bar and a hat, as used in other papers [82,83]. 
In the four-parton antenna functions the hard radiators are "on the edges" and the ure- 
solved particles are "in the middle" . 

Following the general factorisation formula given in eq.(l.l), while omitting the renor- 
malisation and factorisation scale dependences, the contribution of the partonic process 
qq — )• QQ to the total leading order cross section for ti production in hadronic collisions is 

daio = / ^^/.(6)/,-(6)da,,-_QQ, (5.1) 
where the leading order partonic differential cross section is given by 

d(T^q-_^QQ =A^Lod$2(m,P2;P3,P4)|A^4(lQ>2Q,3g-,4g)|2jf^(pi,P2). (5.2) 
The normalisation factor Mlo reads 

where s is the hadronic center of mass energy, (A'^^ — 1) comes from the colour sum, the factor 
l/4Ar^ accounts for the averaging over the spin and colour of the incoming qq pair, and 
(l/2s) is the hadron-hadron flux factor. The coupling is defined as usual: ag = 9'^/47r, C(e) 
is given in eq.(2.27), and C{e) = (47r)'^e~'^'''. This way of expressing the coupling factors, 
with each power of accompanied by a power of C(e), keeps the coupling dimensionless 
in dimensional regularisation. 
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The colour and coupling-stripped amplitude A^4(1q, 2q, 3q, 4^) is related to the full 
amplitude M2(1q, 2q, 3q, 4q) through 

M^ilQ,2Q,3q,Aq) = / (^5i,,iJi,,i, - ^Jn,i,Si3,i?j Ml{lQ,2Q,3q,lg). (5.4) 

(2) 

The jet function denoted by J2 ensures that the top and the antitop are in two separate 
jets. 

The contribution of purely fermionic C(af ) processes to the double real radiation cross 
section for massive quark pair production in hadronic collisions is 

= y" [/«te)/5-(6)d<Tgg-^QQg5- + /g(a)/g(6)dagg^QQ<jg 

+/g(6)/q(6)do-^^_>QQqq] 

+ fqiQfq'{^2)dagg,^QQgg, + ^ ( ) ( ^ ) d^^^/ ^ Q Qg^/ ] |, (5.5) 

with the partonic cross sections given by 

d^qq^QQq'q' = d$4(A;Q, Aiq, kg', % ; P?, P?) |M°-^qq^, | ^ ^ (A;q, kQ, kg>, kq>) (5.6) 

d^gg^QQgg = d$4(fco, kg, kg, kg] Pq, Pg)\M^^-_^QQ^.\'^ J^^\kQ , kg, kg, kg) (5.7) 

dagg,^QQgg, = d 4 ( Zcq , , , ^ ; , ) | M^^ , ^ ^ ^ ^ ^ , | ^ ^ ( , , , /c^ ^ (5.8) 

d^gg^QQgg = d^ ^{kQ , k Q , k g , k g, P g , P g)\M^^^_^QQ^^\^ J^^\kQ , k Q , k g , k g) . (5.9) 

Charge conjugation relates all other partonic cross sections in eq.(5.5) to one of those given 
in eqs. (5.6-5.9). For example dagg_^QQgg is identical to dagg_^QQgg, etc. 

The jet functions appearing in eqs. (5.6-5.9) are all of the type Jg^^ and correspond to 
the selection criteria of a 2-jet event: out of four partons, from which two are a QQ pair, 
an event with two jets is built. Each of these two jets has the heavy quark Q or the heavy 
antiquark Q in it. 

5.1.2 Colour structure 

In order to obtain the subtraction terms, the colour decomposition of the real matrix ele- 
ments present in the partonic cross sections in eqs. (5.6-5.9) has to be performed. However, 
since colour ordering does not distinguish between initial and final state partons, it is suffi- 
cient to consider the colour decomposition for the unphysical processes QQqqq'q' and 
obtain the matrix elements needed for ti production by appropriate crossings. We shall 
follow this strategy below. 

The colour decomposition of the amplitude for — > QQqqq'q' has been derived in [89]. 
For completeness we recall it here. Since there are no gluons, the colour structure of the 
amplitude is rather simple: no adjoint representation indices are present, i.e. no T" matrices 
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appear, and all colour factors simply consist of products of fundamental representation 
deltas. Taking all particles as outgoing, we can write the amplitude for the fictitious 
process QQqqq'q' as 



1 



+ 



, (5.10) 



The sublcading colour amplitudes in eq.(5.10) always contain a colour chain made of a 
quark-antiquark pair of equal flavour. This fact manifests itself in the corresponding colour 
coefficients, which have a factor Si^j^ (with i and j of the same flavour). The presence of 
this factor is related to the fact that this qq pair splits from a photon-like propagator and as 
we shall shortly see, this has important consequences in the way in which the colour-ordered 
sub-amplitudes factorise in their double soft limits. 

For those processes involving two pairs of identical flavour quark-antiquark pairs, the 
colour decomposition can be obtained from eq.(5.10) by using 

M0(1q, 2q, 3g, 4g, 5g, 6g) = M^{1q, 2^,3,, Ag, 5g,,6g,) " M^{1q, 2q, 3g, 6g, 5g, Aq') , (S-ll) 

where the relative minus sign is due to the fact that both amplitudes differ by the inter- 
change of p4 and pq. 

The amplitudes squared for all partonic processes under consideration can be written in 
terms of colour-ordered amplitudes by squaring either eq.(5.10) or eq.(5.11) and crossing the 
corresponding (anti) quarks to the initial state. When squaring, interferences of different 
colour-ordered amplitudes appear. Most of the time, they can, however, be eliminated by 
using the following identities 

M'iilQ,2Q;;3g,4q,;5g',6g') 
= MtilgAq-, ; 3g, 6,-,; ; 5g., 2q) + 7W0(1q, 6,-,; ; 3„ 2^; ; 5,,, %) 

= I [MtilQ, ; 3g, 2^; ; 5,., 6,-0 + M^ilg, ■ 3„ ; 5,., 2q) 

+>t°(lQ, 2q; ; 3„ 6,-.; ; 5,,, %)) (5.12) 

which relate the leading, subleading, and subsubleading colour partial amplitudes. 
5.2 Soft limits of subleading colour contributions 

The procedure to subtract the limits in which a quark-antiquark pair becomes soft is rather 
straightforward for those colour-ordered amplitudes where the neighbours of the qq pair in 
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the colour chain are easy to identify. When the soft pair Cq, dq is colour-connected to the 
hard particles a and 5, the corresponding sub-amplitude factorises as 

\Ml{...,a,Cq--dqA-)\ Saam{rna,rni,)\Ml_2{---.aA---)\\ (5-13) 

where the double soft factor Sacdb has been given in eq.(4.21). It is clear that a suitable 
subtraction term is 

X°(a, Cq, dq, 6)|M°_2(-, a, (5-14) 

since 

X^{a,Cq,dq,b) — > Sacdb{ma,mb), (5.15) 

and the momentum mapping is such that the reduced matrix element in eq.(5.14) gives 
l-^n-2(---' ^' --OP when the qq pair becomes soft. The type of four-parton antenna func- 
tion needed depends on the identity of the hard radiators a and b, and the type of phase 
space mapping depends on whether these radiators are in the initial or in the final state. 

The same strategy cannot be followed with those subleading colour amplitudes of the 
form Mni---, ct; ; dq, Cq), whose colour factors contain Si^^i^. In these cases the gluon propa- 
gator from which the soft quark and antiquark split is photon-like, and the hard radiators 
cannot be simply identified as the pair's adjacent particles in the colour chain, because the 
concept does not even apply. 

In the next section we shall explain how we can build antenna subtraction terms for these 
subleading colour amplitudes by recalling the factorisation properties of colour-ordered 
amplitudes in their soft limits. Since, from a physical point of view, the emission of a soft 
qq pair is closely related to the emission of a single soft gluon, we shall start by addressing 
the problem for this latter case. Note that our results will be obtained in the most general 
massive case, from which the massless counterparts can be inferred by setting all masses 
to zero. 

5.2.1 Emission of a single photon-like soft gluon 

When a soft gluon is emitted between the hard particles a and b, a colour-ordered amplitude 
factorises as 

M%..,a,g,b, ...)| e,ipq,X) [J^{pg) - J^P,)] A1°_i(..., a, 6, ...) (5.16) 
where e^(p5,A) is the soft gluon's polarisation vector, and the soft currents are given by 

V • Pg 

After squaring eq.(5.16) we obtain the well-known formula 

\Ml{..., a, g, b, ...)|2 Sagbima, mb)\Ml_i{..., a, b, ...)|2, (5.18) 

with 

Sagb{ma,mb) = e^(pp, A)e^(pp, -A)^ [JaiPg) - JbiPa)] Kbs) - ^biPg)] 
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(5.19) 



By evaluating the second line in eq.(5.19) explicitly with the current given in eq.(5.17), 
the form in eq.(4.11) for the soft eikonal factor given is obtained. To get this factor, we 
furthermore use the fact that colour-ordered amplitudes are gauge invariant which ahows 
us to replace 

XI ^i^'^Ps' ^)^^{P9^ -dixv (5.20) 

\=± 

in eq.(5.19). 

It should be stressed once more that in this case, the hard radiators a and h are the 
neighbouring particles of gluon g in the colour chain, and they can be both fermions or 
gluons. This statement is no longer true for sublcading colour amplitudes in which the soft 
gluon g is photon-like. In such cases, the abelian nature of the soft particle implies, on one 
hand, that it does not couple to the hard gluons which the amplitude might have, and it 
also implies, on the other hand, that all the fermions involved in the process can act as 
hard radiators. Therefore, the factorisation of a subleading colour amplitude in the limit 
where a photon-like gluon becomes soft reads 

Ml{...,a--gy^ E '^i^Pa)' E ^^(P.) J a), (5-21) 

where {q} stands for all final state quarks and initial state antiquarks in the process, and 
{q} stands for all final state antiquarks and initial state quarks. Squaring eq.(5.21) and 
rearranging the result into a more convenient form, we obtain the following factorisation 
for the squared subleading colour amplitude: 

\Mr,{...,a;-g^)\ — > \ 2^ Sigj{mi,mj) - - Sigj{mi,mj) 

^ ie{(;} (*J)e{9} 
je{g} i^i 

-\ E Sigj{mi,mM\Ml_^{...,a)\\ (5.22) 
(ij)e{g-} / 

This shows that, unlike the case presented in eq.(5.18), the soft factor corresponding to a 
soft photon-like gluon is a combination of several of soft eikonal factors given in eq.(4.11). 

It can be immediately seen from eq.(5.22) that in order to subtract the soft limits of 
parton 5 in a squared sub-amplitude of the form |A^^(..., a;;g)\^^ we can use the following 
combination of antenna functions and reduced matrix elements ^: 

E Xl{i,g,j)\Ml_^{...ra\i,^)\^ 

je{5} 

^Note that this result was already obtained in [89] where it was formulated slightly differently and 
concerned the soft behaviour of interference terms. 
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2 Xlii,9,j)\M' 







l( 




2 ^ Xl{z,g,j)\M. 







n 



l( 



(5.23) 



(»>j)e{g} 



The subscripts in the reduced matrix elements mean that the phase space mapping 
for each term is different, since the hard radiators in each term are different. Since in the 
hmit — 7- the antenna functions X^{i,g,j) give the soft eikonal factors Sigj(mi,mj) 
and the phase space mappings reduce to the Born kinematics, it is clear that eq.(5.23) 
constitutes a suitable subtraction term for eq.(5.22). 

Two crucial properties of the subtraction term in eq.(5.23) should be pointed out. The 
first of these is that it subtracts not only the soft limits of the photon-like gluon, but also 
all the single collinear limits in which this gluon participates. The second property is that 
it does not introduce any spurious singularities that are not present in the real radiation 
matrix elements. 

5.2.2 Emission of a soft qq pair splitting from a photon-like gluon propagator 

When a soft quark-antiquark pair is emitted between partons a and b in the colour chain, 
the colour sub-amplitudes factorise as 



Mn{-,a,Cq;;dq,b,...)\ -U [usa{Pd)lixVsciPc)] {JaiPcPd) - JbiPcPd)) Mn-2{-,a,b,...), 



where the soft currents are given by 

Squaring eq.(5.24) and summing over the spins of the soft quark and antiquark gives 



with the double soft factor Sacdbi'ma, mb) encountered in eq.(4.21). In terms of the currents 
given above, it can be rewritten as 

Sacdbirria, rub) = i^{i)dltj.t>cyv) (JaiPcPd) - JbiPcPd)) (JaiPcPd) - JbiPcPd))^ (5-27) 

By explicitly evaluating the trace in eq.(5.27) and using the definition of the currents in 
eq.(5.25) the expresssion given in eq.(4.21) is obtained. 

When a qq pair becomes soft in a sub-amplitude whose colour factor contains Si^^i^, the 
factorisation given in eq.(5.24) docs not hold. This is due to the fact that the presence of 
Si^^i^ factor indicates that the gluon propagator from which the soft pair splits is photon- 
like. As discussed in the previous section, in this case the soft pair can be emitted from all 



(5.24) 



\M„{...,a,Cq;;dq,b,...)\ 




Sacdbirua, m6)|A1°_2(..., a, 6, ...)p 



(5.26) 
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the other fermions present in the process. Ah of these can act as hard radiators. Thus, in 
the Hmit PcPd these type of sub-amphtudes factorise as 

M%..,a;;dq,Cq)'"'^ [usa{Pd)lnVp^{Pc)]l ^ JtiPcPd)-^ J^{pc,Pd)]M'i_2{-,a). 

V ie{q} je{q} / 

(5.28) 

Squaring eq.(5.28) we find 

\Mn{...,a;;dq,Cq)\'^'' ^ i ^ Sicdj{mi,mj) - ^ ^ Sicdjimi,mj) 



I E S,,djimi,mj)]\M'i_,{...,a)\\ (5.29) 



2 

(i,i)e{g} 

which suggests that a suitable subtraction term for these Umits can be symbolically written 
as 

J2 i^Kh c, d,j) - Xlii, c, d)Xl{ic, cd,j) - Xlic d,j)Xl{i, cd, dj)\Ml_2i-, 
ie{q} 

-I E iX4ihc,dJ) - Xl{i,c,d)Xlirc,cdJ) - Xlic,dJ)Xlii,cd,m^^^ 
-I J2 (^4 c, d,j) - Xiii, c, d)Xl{ic, cd,j) - Xiic d,j)Xlii, cd, ^•)|A^n-2(-, 

ihJMq} 

(5.30) 

where the products of the three-parton antennae X^ remove the single unresolved be- 
haviours present in the four-parton antennae X^. CiTicially, we have checked that this 
combination of antenna functions not only explicitly subtracts double soft limits in sub- 
amplitudes of the form |A^J^(...,a; ;dq,Cq)\'^, but also reproduces the triple collinear limits 
that involve dq, Cq in these sub-amplitudes without introducing any spurious singularities. 

5.3 Fermionic processes for QQ production at LHC 

In the following, we shall explicitly write the real radiation matrix elements for each process 
in terms of colour-ordered amplitudes squared, and give the corresponding double real 
subtraction terms capturing single and double unresolved radiation features. The colour 
decompositions as well as the subtraction terms of single unresolved limits have already 
been derived in [89]. The double unresolved parts are new. 

5.3.1 The qq — )• QQq'q' process 

By squaring eq.(5.10), using eq.(5.12) to remove interferences, and crossing 3g and 4^ to 
the initial state, we can write the amplitude squared for qq — > QQq'q' as 

|M0(1q, 2q, Iq, bq>,Qq>)\^ = /(iV^ - 1) 
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NJ IMU^Q, 4g- ; 3,-, 6gr, ; 5g,,2Q)\' + IMU^Q, 6gr, ; 3g, 2q; ; 5,,, 4,) 



+— ( |A^^(1q, 4,; ; 3,, 2^; ; 5^., 6,-01' + \M"ei^Q, G^r, ■ 3,-, 4,; ; 5,., 2g)|^ 



+|7W1|(1q, 2q; ; 3,-, 6,-,; ; 5q>,4q)f - 3|7W^(1q, 2q; ; 4,; ; 5,,, 6^ 



(5.31) 



This amplitude is singular in the 5g'||6q' single coUinear limit, as well as in the 5q',6q' 
double soft limit and the 4y||5y'||6y' and 3y||5q'||6y' triple coUincar limits. To subtract these 
divergencies from the differential partonic cross section we use the following subtraction 
term 

a.\2C'(e)2 



2ttJ C(e)2 

\eI{Iq, 5,,, 6,-0 {\Ml{{l^)QAq; ; 3,-, (56),, 2q)|2 

+\Ml{{lb)Q, (5~6)5,4,; ; 3,, 2q)\^^ A^\p{^,P2.P^q) 

1 



+ ;^^3(2q, 5,', 6,-0 |A^^(1q, 4,; ; 3,-, (56),, (25)q)|^ 



+|A^^(1q,(56)„4,;;3,-, {2b)Q)Aj^i\pr,nr.,P^%) 



+ ( s2(1q, 6,-,, 5,., 4,) - -El{lQ, 5,., 6,-0^^((15)q, (56),, 4,) 



-^£;3"(2Q,5,.,6,-0A^(lQ,(56)„4,))|>l^((156)Q,2Q,3,-,4,)|2jf (pIi6,P2) 



+ i?^(2Q,5,.,6,-.,3,-) - -i?3^(lQ, 5,., 6,-0^5^(2^,(56),, 3,-) 



-^£;3°(2q,5,,,6,-0^^((25)q,(56)„3,-))|M^(1q,(256)q,3,-,4,)|^J^^^(pi,P^6) 



0/ 



1 2 t(2). 



+ 



Nr. 



^EiilQ, 5,,, 6,-0 (^|A^^((15)q, 2q; ; 3,-, (56),, 4,) 
+|Atg((r5)Q,(5~6)„2Q;;3,-,4,)p 

0^/^K^„ 0_ Q A (Kh\ M2 \ t(3)/ 



-2|A^^((15)Q,2Q,3,-,4„(56)^)njr(;^5,P2,pr6) 



+ 2^3(2q, 5,,, 6,-0 (^|A1^(lQ, (25)q; ; 3,-, (5~6)„ 4,)!^ 

+|A^0(1q,(51))„(2~5)q;;3,-,4,)|2 

-2|7W0(1q, (25)q,3,-,4„ (5l))^)|2') jf )(pi,p^5,pr6) 



54(1q, 6,-., 5,., 2q) - -El{lQ, 5,,, 6,-0^^((15)q, (56),, 2q) 

-^E3°(2q,5,,,6,-0^^(1q, (5~6)„ {2h)Q)\ |M2((156)q, (256)q, 3,-, 4,)|Vf (pITe.P^e) 
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- (^B2iSg, 6g,,5g>,4g) - ^ElilQ, 5,,, 6,-0^^(%, (56),, 4,) 
-2(^s0(1q,6,-,,5,.,4,) - ^ElilQ,5q,,6q,)Aliir5)Q,{56)g,lg) 

-2(^Bi{2Q,5g,,6g,,3g) - ^El{lQ,5g,,6g,)A^s{2Q,{m)g,3g) 

-^El{2Q,5g>,6g>)Al{i25)Q,i56)g,3g)yMl^^^ 
+2 (^B^ilQ, 6,-,, 5,,, 3,-) - ^ElilQ,5g,,6g,)Ali{15)Q, (56),, 3,-) 

-IeI{2q, 5g>,6g>)Al{lQ, (56)g, 3,-) j |7W2((156)q, 2q, 4,)|2 (^1^6,^2) 

+2 (^S°(2q, 5,., 6,-., 4,) - ^ElilQ, 5,., 6,-0^§(2q, (56),, 4,) 

-^^3^2Q,5,^6,-0^^((2~5)Q,(5~6),,4,))|>tO(lQ,(256)Q,3,-,^ |, 

(5.32) 

where the pieces containing the ^4 antennae subtract the double soft and triple coUinear 
limits, while those multiplied by E^ antenna functions subtract the single collinear limit 

5g/||6g'. 

To subtract the double unresolved limits of the subleading colour sub-amplitudes squared 

|7Wg(lQ,4g;;3g,2Q;;5g/,6y/)p and |7W|](1q, 2^; ; 3^, 4^; ; 5,/, 65/)^ we followed the proce- 
dure described in section 5.2.2 and used eq.(5.30). The generic four-parton antennae 
appearing in cq.(5.30) are here given by type antennae, while the generic X^X^, also 
present in eq.(5.30), are given above by products of £"3 and Ag-type antennae. 

5.3.2 The qq QQqq process: identical quark flavour contributions 

The amplitude squared for the identical flavour process qq QQqq can be written accord- 
ing to eq.(5.11) as 

IMe^lQ, 2q, % 4,, 5,, 6,-)p = |MO(1q, 2q, %, 4,, 5,^ + |M60(1q, 2q, % 

+2Re (M°(1q,2q,%,4„5,,,6,-0MO(1q,2q,3,-, V,5,.,6,-)t) . (5.33) 

The Feynman diagrams contributing to this process can be divided into two categories: on 
one hand we have diagrams that are also present in the non-identical flavour case, which 
are included in the amplitude Mg (1q, 2q, 3y, 4g, 5^', 6^/), and on the other hand wc have 
diagrams in which the initial state quark and antiquark go through to the final state. These 
latter kinds of diagrams are specific to the identical flavour case and are included in the am- 
plitude M^(1q,2q,%, 4g/,5g/,6g). The amplitude squared |M^(1q, 2q, 3g, 4q, 5g, 6g)p has 
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therefore a piece that is the same as the non-identical flavour contribution (|Mg (1q, 2q, 3,, 
4g, 5,/, 6g')p), and a remainder piece which we call "identical- flavour-only" . The singular- 
ities of the former have been already subtracted with (5.32) when we also included the 
identical quark flavour possibility in the sum over flavours, so we shall now focus on the 
identical-flavour-only piece 

|MO(1q,2q, 3,-, 4„5„ 65)1^0 = |M0(1q,2q,3,-, V,5,s6,-)|2 

+2Re (mO(1q, 2q, %, 4g, ,6g,)M^{lQ, 2q, 3,-, V, 5,,, 6,-)t) • (5-34) 

In terms of colour-ordered amplitudes, this matrix element squared is given by 

\MlilQ,2Q,3gAg,5g,Gg)\h = 9'{N! - 1) 
x!^N,(^\MtilQ, 6,-; : 3,-, V; ; 5,,, 2q)|2 + |>10(1q, V; ; 2^; ; 5g,,6g)f^ 

-2Re(A^[](lQ, 4,.; ; 3,, 2^; ; 5,,, 6,-)A^[!(1q, 4,; ; 3,-, 2^; ; 5,., G,-/)^) 
-2Re(7W0(lQ, V; ; 6,-; ; 5,,, 2q)A^0(1q, 4,; ; 6^-,; ; 5,,, 2Q)t) 
-2Re(A^0(lQ, 6,-; ; 3,, 2^; ; 5,,, 4,0A^6(1q, 6,-.; ; 2^; ; 5,., 4,)t) 



-2Re(A^0(lQ, 6,-; ; 3,-, 4,,; ; 5,,, 2q)>J0(1q, 65, 
+2Re(A^0(lQ, 2q, ; 3,-, 6,-; ; 5,,, 4,0A^!](lo, 2q 



;%,4,;;5,,,2Q)t) 
;3Q,4g;;5,/,6g/)^) 



+ ^ (^\MI{1q, 6,-; ; %, 2^; ; 5,., 4,0|' + |Mg(lQ, V; ; % %; ; V> 2q)|2 
+ |A10(1q, 2q; ; 4,,; ; 5,,, 6,-)|2 - 3|A10(1q, 2q; ; 3,-, 6,-; ; 5,,, V)P 

+ ^ (^6Rc(A^0(1q, 2q; ; 3,-, 6,-; ; 5,., V)X^(1q, 2q; ; %, 4,; ; 5,,, 6,0^) 
-2Re(A1^(lQ, 2q; ; 3,-, 6,-; ; 5g, ,4g,)MtilQ, 2q- ; 3,-, 6,-.; ; 5,., 4,)t) 
-2Re(M°(lQ, 2q; ; 3,, 4,,; ; 5,., 6,-)A^^(1q, 2q; ; %, 4,; ; 5,,, 6,-.)^)) } . (5.35) 

The amplitude squared for the identical-flavour-only contribution is singular in the 
initial-final collinear limits 4g'||5g' and 3,1 |6g, in the double collinear limit 4g/||5g' -|- 3,1 |6g, 
as well as in the triple collinear limits 3^||5,/||6g and 4,/||5,/||6g. The interference terms in 
eq.(5.35), which are interferences between partial amplitudes coming from Mg(lQ,2Q,3^, 
4<j,5q',6^') and Mg (1q, 2q, 3,, 4,/, 5,', 6,), are only divergent in the triple collinear limits. 
The full subtraction term for this contribution is 



^^qq^QQqq,lO = ^LO ) -^d^^ipi, P2, P5, PG', P3, Pi) 



X<!iVe 



\Ei{lQ, 5,,, V) (|>1^((15)q, I„ 6,-; ; 3,-, 2q)\' 



+|A^^((15)q, 6,-; ; 3,-, 4^, 2q)\' j ipr5,P2,P6) 
+IeI{2q, 5,,, 4,0 (\M'i{lQ, tg, 6,-; ; 3,-, (2~5)q)|2 
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2\ 7(3), 



+-EI{2q, 6,-, %) 3„ 4,,; ; 5,,, (26)^)^ 



+|>l0(lQ,4,,;;5,.,3g,(26)Q)|2jjf (pi,p^6,P5) 
+ ('50(65,5g,,4,,,3g-) - |£;3°(lQ,55.,4,0^^(6g-,lg,3g-) 



+ (^0 (5,., 6,-, 3,-, 4,0 - ^£;0(lQ,6g-,3,-)^0(5,,,3„4,0 



--^3°(2q, 6,-, 3,0^°(5,^ 3g, 4,0 ) |A^^(iQ, % J,-^ 



--ElilQ, 5,,,ig')Eli2Q, 6,-, %) (^|>t^((15)Q, tg, (26)q)|2 
+|7W2((r5)Q,l3,3„(2l))Q)rtjf)(pl5,p56) 



--£;3°(1q, 6,-, 3g)Ei{2Q, 5,,, 4,0 |A^2((16)q, 3^, ig, (25)q 



+ 



+|A^°((16)Q,4„3g,(25)Q)pjjf (pl6,P~25) 

+2C0(3,-, 5„ 65, 4,)|>jO(iQ, 2q, §5, \)\^4^\pi,p~2) 
+2C2{lg, 6g, 5„ 3,-)|>l2(iQ, 2q, K)\^4^\p~i,P~2) 
ElilQ, 5,,, 4,0 f |Xg((r5)Q, 4„ 2q; ; 3,-, 6,-)|' 



+ |A^°((15)q,2q;;3,-,I„6^ 



-2|7Wg((15)Q,2Q,3,-,6,-,l^)|2)jf (pr5,P2,P6) 



+-£;30(2q, 5,,, 4,0 (^|A^g(lQ, t„ (25)q; ; 3,-, 6,-)p 
+|>tg(lQ,(2~5)Q;;3,-,4„6,-)P 
-2\Mt{lQ, (2~5)q, 3,-, 6,-, 4^)p') 4^\pi,m5,P6) 



+-ElilQ, 6,-, 3,-) (^|M0((16)q, 3^, 2q; ; 5,,, 4,, 

+|A1°((r6)Q,2Q;;5,.,3g,4,0|=^ 
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-2|A^^((16)q, 2q, 5,., V, 3^)|2 ) (pr6,P2,P5) 



+ 2^3 (2q, 6,-, 3,-) 3„ (26)q; ; 5,,, V)r 

+|>t«(lQ,(2~6)Q;;5,,, 3^,4,01' 

s2(6,-,5,.,4,,,3,-) - ^i?3°(lQ,5,,,4,0A°(6,-,4„%) 
1 



-^£;3^2q, 5,^ 4,0^^(6,-, 4^, 3,-) ) 2q, 3,-, 4,)| (pi,p~2) 



- { 5^(5,,, 6g-, 3g-, 4,0 - -^3^(1q, 6g-, 3g-)Al^(5,,, 3g, 4,, 



-^£;3^2q, 6,-, 3,0A^(5,^ 3^, 4,0 ) |A1^(1q, 2q, 3,-^ 4,01' (pl,P2) 



+^£;3"(1q, 5,^ 4,0£^3"(2q, 65, 3,0|>1^((15)q, 3^, 4^, (26)q)|2 )(pl5,p56) 



+^E^{1q, 6g, h)Ei{2Q, 5,., V)|M^((r6)Q, 3^, 4^, (25)q)| (pl6,p~25) 



iV2 



202(3,-, 5„ 6,-, lq)\M\{lQ, ~2q, %, \)\^J?\pi,P2) 
+2Cl{lq, 6,-, 5„ 3,0|A^5(iQ, 2q, 3,-, )(pi,:p2) 



(5.36) 



The choice of one of the hard radiators in the 5° antennae used to subtract the triple 
collinear Hmits is rather ambiguous. Any of the three particles that are not involved in 
the limit can be chosen as a radiator, since the only role it plays is to recoil and ensure 
overall momentum conservation. We decided to use always antennae with massless hard 
radiators. 

5.3.3 The qc^ — ^ QQqq' process 

By squaring eq.(5.10), using eq.(5.12) to eliminate the interferences terms, and crossing 4, 
and 6q' to the initial state, the amplitude squared for the process qq' — >■ QQqq' in terms of 
partial amplitudes squared is obtained. It reads 

|MO(1q, 2q, 3„ 4„ 5,., 6,01' = g\N^ - 1) 

Nc(^\Ml{lQ, 4,; ; 3„ 6,,; ; 5,,, 2q)|2 + \MI{Iq, 6,,; ; 3„ 2q; ; 5,,, 4,)|2 

+ ^ (\Mt{lQ, 4,; ; 3„ 2q, ; 5,., 6,0|' + |A^^(1q, V; ; 3„ 4,; ; 5,., 2q)|2 

+\MI{1q, 2q, ■ 3„ 6,.; ; 5,,, 4,)|2 - 3|M°(1q, 2q, ; 3„ 4,; ; 5,., 6,0|') | • (5.37) 

The singularities of this amplitude are: the 3,||4, and 5,'||6,' collinear limits, the 3,||4, + 
5,/||6,' double collinear limit, and the 3,||4,||5,' and 3,||5,'||6,/ triple collinear limits. To 



-45 - 



subtract these divergencies from the partonic differential cross section we use the following 
subtraction term 



,=Mlo[ 



a,. 



^d$4(pi,^'2,P3,P5;P4,P6) 



27r/ C(e)2 

\eI{Iq, 3„ 4,) {\Ml{{Vi)Q, l„ %r, ■ 5,,, 2q)|2 

+|A^°((r3)Q,6,,;;5,.,4„2Q)|2')jf (pr3,P2,P5) 

1 



+;^^3 (2q, 3,-, 4,-) 4„ 6,,; ; 5,,, (23)q)|^ 



+|A^^(1q, 6g.; ; 5,,, 4^, (23)q)|2 ) 2^3,^5) 



+-E^^{lQ, 5,,, 6,0 |7W^((15)q, 6„ 4,; ; 3„ 2q)|^ 



+|>t^((15)Q, 4,; ; 3,, 6^, 2q)\^ ] jf (pr5,^>2,P3) 



+ 2^3 (2q, 5,., 6,0 (^|A^^(1q, 6„ 4,; ; 3„ (25)g)|^ 

+|A10(1q, 4,; ; 3„ B,, (2'5)q)|2') (^,^^5,^3) 



+ ( S^(3„ 6,., 5,., 4,) - -£;3"(1q, 5,., 6,0^^(3,, 6^, 4,) 

-^£;30(2q, 5,., V)^^(3„ 65, iS) \MI{Iq, 2q, 6,-, I,)! (pi,p2) 



+ { i?0(5,,,4„3g, V) - 2£;30(lQ,3g,4,)A0(5,,,4^,6g' 



-^£;3"(2^,3,,4,)A^(5,^45,6,0)|>t^(lQ,2Q,%^6,0|Vf (pi,P2) 



-;^^3 (iQ, 3„ 4,)£;3^(2q, 5,,, 6,0 |A^^((13)q, 4„ 6„ (25)q)| 



+|A^^((13)q,6„4„(25)q)|2 (pl3,p55) 



■-£:30(1q, 5,,, 6,0^3°(2q, 3„ 4,) ( |>t^((15)Q, 4„ 6„ (23)q)|2 



+|>l^((15)Q,6„4„(23)Q)|2jjf (pl5,p53) 

£;3°(1q, 3„ 4,) {\Ml{{Vi)QX^ 2q; ; 5,,, 

+\Mlm)Q,2Q--bq,,lgA')\^ 
-2|X0((r3)Q,2Q,5,.,6,,,4^)pV?^(^6,P2,P5) 



+ 
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+IeI{2q, 3„ 4,) [iMlilQ, ig, (23)q; ; 5,,, 6,01' 

+|A^0(1q,(23)q;; 5,., 23,6,01' 
-2\Ml{lQ,{23)Q,5g>,6g,,t^)\^)4^\pi,m3,P5) 

2 



+ -£;3^(1q, 5,. , 6,0 ^|A^^((15)q, 6g, 2q, ■ 3„ 4, 

■(0/'/'i1r\_ 0_..Q « /I M2 



+|A1^((15)q,2q;;3„6^,4,-) 
-2|M°((r5)Q,2Q,3„4„6^)P^jf (pr5,P2,P3) 

+^^3 (2q, 5,,, 6,0 63, (2~5)q; ; 3„ 4,)|2 

+|A^°(1q,(2'5)q;;3„65,4,)|2 
-2|7W0(1q, (2'5)q, 3„ 4„ 6^)|2 j jf (;,i,p^5,P3) 

-^£;3°(2^,5,,, 6,0^^(3,, 6p,4,)) |>t°(iQ,2Q,6,-,4,)|Vf (pi,P2) 
-(^B^(5,.,4„3„6,0 - lElilQ,3g,4,)Al{5g>,tg,6g>) 

-^Eli2Q,3q,4,)Al{5,,,tg,6q,)^ 2q, I,-,, 6,01' (r,P2) 

+^£;3^1q, 3,, 4,)£;30(2q, 5,^ 6,0|A^2((r3)Q, t^, 6^, (2~5)q)| (i?l3,p~25 



(5.38) 

As for the previous process only massless four-parton antennae whose integrated form is 
available [32,80,81] have been used. 

5.3.4 The qq QQqq process: identical quark flavour contributions 

As we did for the qq — > QQqq process, we split the amplitude corresponding to the qq — > 
QQqq process into a piece that contains the same diagrams as the non-identical quark 
flavour case, and another piece which consists of diagrams which are specific to the identical 
flavour case. After squaring eq.(5.11) we obtain 

rO/-i_ r)_ Q A K cm2 ^ f I ]\/rO t 



|M6^(lQ,2Q,3,,4,,5,,6,)|^ = -(^|M6^(lQ,2Q,3,,4,,5,^6,0r + |M6^(lQ,2Q,3,,4,,,5g,,6g)|^ 
-2Re (mO(1q, 2q, 3„ 4„ 5,., 6,0)mO(1q, 2q, 3„ 4,,, 5,,, 6,)t) ^ . (5.39) 
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The factor 1/2 is a symmetry factor that accounts for the presence of identical particles in 
the final state. The interference term comes with a minus sign since the amplitudes differ 
by the interchange of the initial state quark labels. The identical-flavour-only amplitude 
squared can be written in terms of colour-ordered amplitudes as: 

|MO(1q,2q,3„4„5„6,)|2o = g\Nl - 1)^ 
X |iV, {\Ml{lQ, 6,; ; 3„ V; ; 5^,, 2q)|2 + |A^°(1q, V; ; 3^, 2^; ; 5^,, 6,)!^ 

+2Rc(M[!(1q, V; ; 3^, 2^; ; 5^,, 6g)A^[!(lQ, 4^; ; 3„ 2^; ; 5,,, 6,0^) 
+2Ro{MI{Iq, l^r, ; 3„ 6g; ; 5^., 2q)7W[!(1q, 4,; ; 3„ 6,.; ; 5,,, 2Q)t) 
+2Re(A^0(lQ, 6,; ; 3^, 2^; ; 5,,, V)A^[!(1q, V; 5 3„ 2^; ; 5,,, 4,)t) 
+2Re(7W[](lQ, %■ ; 3„ 4,,; ; 5,., 2q)A^0(1q, 6,.; ; 3„ 4,; ; 5,,, 2q)'^) 
-2Re(7Wg(lQ, 2q; ; 3^, 6^; ; ,%,)MI{Iq, 2q- ; 3g, 4^; ; 5,., 6,0^) 

+ ^(|A^6(1q>6,;;3,,2q;;5,,,4,0P + |K(1q>V;;3.,6,;;5,,,2q)|2 

+|A^0(1q, 2q; ; 3^, 4^,; ; 5,., 6g)|2 - 3|7Wg(lQ, 2^; ; 3„ 6,; ; 5,., 4^01" 

(^6Rc(A^0(1q, 2q; ; 3„ 6,; ; 5,., V)A^1](1q, 2q; ; \, 4,; ; 5^,, V)^) 
-2Re(A^0(lQ, 2q; ; 3,, 6,; ; b^, ,%,)MI{Iq, 2q- ; 3^, %r, ; 5,,, 4,)t) 
-2Re(>t°(lQ, 2q; ; 3,, 4,,; ; 5,,, 6,)AtO(lQ, 2^; ; 3„ 4,; ; 5,,, 6,0^)) \ ■ (5.40) 



Apart from the interference terms, this amplitude squared is the same as that of the non 
identical flavour quark case with p4 and pq interchanged. The interferences are singular in 
the 3g||4g||5g and 35II55II65 triple coUinear limits, and these singularities are subtracted with 
C-type four-parton antennae. Thus, we can write our subtraction term for the identical- 
flavour-only case as 

d%-^QQqq,IO = ^^(^qq'^QQqq'iPi ^ Pg) " 3™ ) ^7^d^i{Pl, P2, PS, P5; P4, Pg) 



1 



C{ef 

<on„ o . «_ ^ m2 7(2),' - 



1 - \ 2Cl{Aq, 6q, 3q, 5,)|A^^(1q, 2q, 6,-, Aq)\'4 'iP~l,P~2) 

+2C2{6q, 4„ 3„ 5,)|M2(iQ, 2q, tq, )(pi,p2) 



(5.41) 

5.4 Infrared structure of the subtraction terms 

In all partonic processes considered above, examining the terms which capture the single 
unresolved limits in the respective subtraction terms (which belong to the subtraction term 
'^^nn\o Si'^cin in eqs. (2.6-2.8)) which are made of three-parton antennae multiplied by five- 
parton reduced matrix elements, one can determine how the infrared poles that arise in 
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these terms upon integration over the appropriate one-particle antenna phase spaces will 
have to cancel against the poles present in the real-virtual fivc-parton contributions. 

Looking at each process individually, we find that only in the subtraction term related 
to the process qq — > QQq'q' the poles generated after integration over the antenna phase 
space will have to cancel against explicit poles present in the one-loop contributions at 
the five-parton level. These poles have to cancel against explicit poles related to the 
process qq — > QQg at one-loop. Indeed, in this case the subtraction term at the six-parton 
level involves the same process (i.e. qq — )■ QQg at tree-level) multiplied by a massive E- 
type three-parton antenna. This antenna captures the final-final quark-antiquark coUinear 
features and its integrated form is related to the massive 

1(1) 

-type operator Iq^ defined 

in [89] and quoted in eq.(A.16). These statements are valid both for the leading and 
subleading colour contributions. 

For all other processes, the e poles in the integrated dcr^'^^^ subtraction terms will have 
to cancel against poles present in mass-factorisation countertcrms. This can be seen from 
the fact that the da^'^^Q pieces of those subtraction terms only contain initial-final E-type 
three-parton antennae. These capture the massless initial-final single collinear features of 
the real radiation matrix elements squared and the only e pole present in their integrated 
form is proportional to the splitting kernel p° , as it can be seen from eq.(A.18). 



6. Results 

For all processes considered in section 5 we have tested how well the subtraction terms ap- 
proximate the double real matrix elements in all their single and double unresolved regions 
of phase space, so that their difference can be numerically evaluated in four dimensions. 
In this section we describe how this testing procedure is performed. 

To verify how well the subtraction terms approximate the double real contribution, 
we have used RAMBO [98] to generate phase space points in the vicinity of the singular 
regions and computed the ratio 



'NNLO 



R = (6.1) 
do-? 



for each of these points. As before, da^^^Q stands for the double real radiation corrections 
while dcj^jy^Q is the corresponding subtraction term. In each unresolved limit we defined 
a control variable x that allowed us to vary the proximity of the phase space points to the 
singularity. For the difference da^^^^o ^ ^^^NNLO finite and numerically intcgrablc in 

four dimensions, the ratio R should approach unity as we get close to any singularity. The 
phase space points were generated with a fixed centre-of-mass energy of y/s = 1000 GeV, 
the heavy fermions were given a mass of 174.3 GeV, and the two hard jets were required 
to have pr > 50 GeV. 

In the following, wc shall present our results for one particular partonic process, namely 
for the qq — )■ QQq'q' process which we separate into identical and non-identical quark 
contributions. 
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6.1 The qq — )• QQq'tf process 

We start by presenting the behaviour of the subtraction terms in double unresolved limits, 
before describing how those approximate the real matrix elements in single unresolved 
regions of phase spaces. 

6.1.1 Double soft limit 

For the process qq — )• QQq'q' , the double soft phase space configurations are characterised 
by the QQ pair taking nearly the full center of mass energy of the event s, leaving the 
massless final state qq pair with almost zero energy as depicted in fig. 1(a). In fig. 1(b) we 



5000 




(a) (b) 

Figure 1: (a) Ilustration of a double soft event, (b) Distribution of R for 10000 double soft phase 
space points. 

show the ratio between the double real radiation matrix element and the subtraction term 
for three different values of a; = (s — si2 — 2mQ)/s. It can be seen that as the soft qq 
pair takes a smaller share of the total energy, i.e. as x becomes smaller, the peak of the 
distribution around ii = 1 is sharper. This is a sign that the approximation improves as 
the limit is approached. 

6.1.2 Triple collinear limits 

The other double unresolved configurations in which the amplitude for this process is 
singular are the triple collinear limits 3g||5q'||6g' and 4g||5q'||6g/. In fig. 2(b) we show how, 
as we make x = —s^^q/s smaller, that is, as we get closer in phase space to the singularity 
of the real radiation matrix element, the histogram becomes more pronouncedly peaked 
around unity, signaling again that the approximation is good. Similar results which are 
not shown are obtained for this process for the triple collinear limit involving parton Aq. 

6.1.3 Single collinear limit 

The only single unresolved limit in which the matrix element for this process is divergent 
is the final-final collinear limit 5q'||6g'. The topology of this kinematical configuration is 
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Figure 2: (a) Topology of the triple coUinear limit (b) Distribution of R for 10000 double triple 
collinear space points. 




(b) 



Figure 3: (a) Kinematics of the final-final single collinear limit, (b) Distribution of R for 10000 
single collinear phase space points. 



shown in fig.3(a), and the divergence is approached as the ratio x = s^e/s becomes closer 
to zero. 

As it was discussed in section 4.4, the azimuthal terms that the antenna subtraction 
term develops in its single collinear limits and those acquired by the real radiation matrix 
element do no match. This is the reason why the distributions in fig. 3(b) are not very 
sharply peaked, and this is also the reason why the peaks do not become sharper as x 
becomes smaller. As explained in section 4.4 these angular terms cancel, both in the matrix 
element squared and in the subtraction term, when a point in the collinear region of phase 
space is combined with another point that differs from the original one by a rotation of 7r/2 
around the collinear axis. In fig. 4 we see how the ratio of the real radiation matrix element 
and the subtraction term is substancially improved when the aforementioned combination 
of phase space points is performed. 
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Figure 4: Distribution of R for 10000 single coUinear phase space points after azimuthal terms 
have been ehminated by combining phase points related by a 7r/2 rotation about the coUinear axis. 

6.2 The qq — )• QQqq process: identical quark flavour contributions 

When the incoming qq pair is of the same flavour as the outgoing one, the corresponding 
amplitude becomes singular in more regions of phase space, and new diagrams which are 
divergent in the regions where the non-identical flavour contributions were not singular are 
introduced. As we did in the previous section, we shall consider the identical-flavour-only 
piece, which contains the triple coUinear limits 3g||5q||6g and 4q||5g||6g, the double coUinear 
limit 3g| |6g + 4g| |5g and the initial- final single coUinear limits 3q||6q and 4q||5g. 

6.2.1 Triple coUinear limits 

The only limits shared by the identical and non-identical quark flavour contributions are the 
triple coUinear limits 3g| |5g| |6g and 4g| |5g| |6g, although the splitting functions corresponding 
to these limits in each case are different. 
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Figure 5: Distribution of R for 10000 triple coUinear phase space points. 
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As it can be seen in fig. 5, the ratio between the real radiation corrections and the sub- 
traction term converges to unity as the limit is approached. Again, similar results which 



are not shown are obtained for the hmit involving the incoming quark 4^. 



6.2.2 Double collinear limit 

Another important feature of the identical-flavour-only contributions is that, containing 
amplitudes in which both incoming massless particles go through to the final state, double 
collinear limits are possible as depicted in fig. 6(a). It can be seen in fig. 6(b) that, in the 





(b) 

Figure 6: (a) Topology of a double collinear limit with two initial final collinearities. (b) Distri- 
bution of R for 10000 double collinear phase space points. 

double collinear limit, the presence of azimuthal terms also spoils the convergence of the 
subtraction term to the real radiation corrections. These terms can be averaged out by 
combining four points in phase space related by two 7r/2 rotations about each of the two 
collinear directions respectively. The improvement of the convergence when these angular 




Figure 7: Distribution of R for 10000 double collinear phase space points with azimuthal terms 
averaged out. 



averaging is carried out can be clearly seen in fig. 7. 



6.2.3 Single collinear limits 

As opposed to the non-identical flavour case, in which there was only one collinear limit 
involving the massless final state qq pair, the identical flavour contributions contain two 
initial-final collinear limits: 3g||6q and 4q||5g. The topology of these unresolved limits 
is illustrated in fig. 8(a). In fig. 8(b) we show how, in the presence of angular terms, the 




(a) (b) 

Figure 8: (a) Topology of an initial-final single collinear event, (b) Distribution of R for 10000 
phase space points in the single collinear region. 

distributions for the ratio R are not peaked around unity. As shown in fig. 9, the convergence 
of the subtraction term is again improved when we combine two events that differ by a 
rotation of 7r/2 of the collinear pair about the collinear axis. This combination of phase 
space points cancels out the azimuthal correlations, as it did in the single collinear limit 
discussed in section (6.1.3). However, it should be mentioned that for initial-final limits 
this procedure is not as straightforward as it is in the final-final case, since the rotation 
of the collinear pair takes the incoming particle involved in the collinear limit out of the 
beam axis ^. 

Similar results were obtained for the 4g||5q collinear limit. 

The same checks have been performed on the subtraction terms (5.38) and (5.41) cor- 
responding to the processes qq' — )• QQqq' and qq — )• QQqq respectively. The convergence 
of these subtraction terms to their corresponding squared matrix elements has been tested 
in the same way as described above, and the numerical results obtained are very similar to 
the ones shown in the histograms of this section and related to the processes qq — )• QQq'q' 
and qq — )• QQqq. 

7. Conclusions 

We presented the extension of the antenna formalism required for the calculation of hadronic 

^To compensate for this, the rotated phase space trajectory has to be boosted to a frame in which the 
incoming particles are back-to-back, and then another rotation is needed in order to bring the incoming 
momenta back to the beam axis. 
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Figure 9: Distribution of R for 10000 single coUinear phase space points after azimuthal terms 
have been ehminated by combining phase points related by a 7r/2 rotation about the coUinear axis. 



processes involving a pair of massive final states in association with jets at the NNLO level. 
We applied this subtraction framework to the computation of the double real radiation cor- 
rections to the hadronic tt differential cross section arising from partonic processes with 
fermions only. In section 2, we outlined the general structure of the subtraction terms which 
capture all single and double unresolved features of the real matrix elements squared, while 
focussing on the differences caused by the presence of massive final states compared to the 
NNLO antenna formalism developed for the computation of massless QCD reactions at 
hadron colliders [82,83]. The genuine unknown NNLO tree-level massive four-parton an- 
tenna functions required to capture the double unresolved singular behaviour of the real 
matrix elements when the two unresolved partons are colour-connnected, are derived and 
their limiting behaviour is presented in sections 3 and 4 respectively. The subtraction 
terms given in leading and subleading colour contributions are constructed in section 5 
for all partonic processes involved. The subtraction of double soft singularities from a 
fermion-antifermion pair in subleading-colour contributions requires a special treatment 
which is also explained in that section. Finally, in section 6 we checked the validity of our 
subtraction terms numerically by verifying that they approximate the real matrix elements 
in all the single and double unresolved configurations in a point-by point manner. The 
numerical results showed that the combination of antenna functions multiplied by reduced 
matrix elements present in the subtraction terms correctly describes the infrared singular- 
ity structure of the real matrix elements in all limits including also single collinear limits. 
This is achieved provided the azimuthal terms associated with these collinear limits are 
correctly handled. 

The double real subtraction terms related to the all fermion processes presented here 
provide a substantial step towards the calculation of the NNLO corrections to top quark 
pair production at the LHC. Future steps include in particular the computation of the 
remaining double real subtraction terms related to partonic channels involving gluons and 
the computation of mixed real-virtual contributions for all partonic channels involved. 



-55- 



8. Acknowledgements 



We would like to thank Joao Pires for many useful and stimulating discussions. This 

research was supported by the Swiss National Science Foundation (SNF) under contract 
PP0022-1 18864 and in part by the European Commission through the 'LHCPhenoNet' 
Initial Training Network PITN-GA-2010-264564', which are hereby acknowledged. 



A. Three-pcirton antennae 

In this appendix we list the three-parton antenna functions used in the subtraction terms 
of section 5 together with their single unresolved limits. The massive antennae that we 
give here can be found in [88,89], while massless ones are given in [32,79]. All the single 
unresolved factors can be found in section 4.1. 

Two types of antennae are needed: E-type antenna functions, which are used to subtract 
single and double collinear limits in all the different pieces of the subtraction terms, and A- 
type antenna functions, which are needed in the diT^^^^ pieces since they are the reduced 
matrix elements to which antennae collapse in their spurious single collinear limits. 

A.l A- type antennae 

The only final-final A- type antenna that we use in our subtraction terms is massive. It has 
been computed and integrated in [88], and its unintegrated form is 

.0/1 o 9 N 1 f 2gf2 , 2S12 2S12 S23 Sl3 

^sK^Q^-^g^^Q) = 7 V ' ' ' ' 

(E^„, + 2mlj V«13S23 S13 S23 Sl3 ^23 



(- 



2 / 8^12 2512 2512 2S23 2 2 ^ 2.Si3 

^ V«13S23 Si3 ^23 ^IZ «13 •S23 §23 



^23 *13 



with E'^^ = {pi +P2 H-ps)^- This antenna has only a single soft limit: 

Al{lQ, Zg, 2q) 5i32(mQ, mq). (A.2) 
The flavour violating initial-final antenna reads 

.0/1 q 9 A 1 ( 2.sf2 , 2S12 2si2 Sl3 S23 
^3(lQ>3g,2q) = ^- "Y — — + — + — + — 



(q2 _ ml) 



\'Sl3S23 Si3 S23 S23 Sl3 



2S12 


2S12 ^ 




Sl3 


S23 


^)) 


+ 0(e), 





+ (A.3) 



'13 ■'13 

where = —{pi — P2 +^3)^- The integrated form of eq.(A.3) can be found in [89]. This 
antenna has a soft gluon limit as well as an initial-final collinear limit: 

A0(1q, 3^, 2g) S^32{mQ, 0), (A.4) 
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Two different initial-initial A-type antenna functions are used in our subtraction terms. 
They are given by 

AHi„3„2,) = ^J-^ + ^-^ + '-^-'^-''^)+ OM. (A.6) 

\ S13S23 Si3 523 S23 SlS/ 
\ S13S23 Si3 S23 S23 513/ 

with = —{p3 — pi — P2)'^ in eq.(A.6) and = —{pi — P2 — Pd)^ in eq.(A.7). These 
antennae have been integrated in [79]. Their infrared limits are 

A°(i,-,3„2,)'-^°5i32, (A.8) 
AUk, 3g, 2,) —P,g^Q{z), (A.9) 

Sl3 

A'iik, 3„ 2,) J-P^^^q(^), (A.IO) 
S23 

and 

Alilg, 3g, 2g) ^P^^^^^z) (A.ll) 

Sl3 

A. 2 E-type antennae 

For the subtraction of final-final coUinear limits we use 

El{lQ,3q, Ag) = ^ -2 f si3 + si4 + ^ + ^ - 2E,mmQ] + 0(e), (A.12) 

(el -ml) V ^34 .34 J 

where = {pi + p3 +^4)^. The integrated form of this antenna can be found in [88], 
and its only infrared limit is 

ElilQ, 3„ 4,-) J-P ^c(z) (A.13) 
S34 

On the other hand, to subtract initial-final collinear singularities we use the following 
initial-final antenna 

£;30(1q,3„4,) = -- i -2 (~su + Si3-'^-'-^-2mQm^)+0{e), (A.14) 

(Q^ + ml) V ^34 SS4 J 

with = — (P1+P3— P4)^, and = \/ —Q^ is the mass of the incoming neutralino denoted 
by X and present in the process x<? — ^ Q(l defining this initial- final £^3 antenna [66]. The 
only infrared limit that this antenna has is 

El{\Q,\Ag) -^Pgg^Qiz). (A.15) 

S34 
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A. 3 Integrated forms of massive antennae 

To supplement the discussion given in section 5.4 on the infrared structure associated with 
subtraction terms capturing single unresolved limits, we list the integrated forms of the 
two three-parton antenna functions that we use in those pieces of the subtraction terms. 

After integration over the corresponding 1 — > 3 antenna phase space, the massive final- 
final antenna function reads [88] 



6(1 - ^2)3 

-2/u'^(-3-3/x + At^)ln(/n2 



6 + 3/x - 1 V + 14/ - 3/i^ - 14/x2 + - Sfi^ - Qfi^ 



+ 0{e) 



where the infrared singularity operator is given by [86, 89] 



[(1) 



,Nj. €,SQg,mQ,0, 



1717 



SQg + '<3 2r(l - e) 



The integrated form of the massive initial-final -Eg antenna is given by [89] 



^3(4,;3„1q) 



r(i - 6) 



+ — (2 -2x + x^) [-2 + 21n(l - x) - ln(x) - ln(l - xqx)\ 



3a; - 2x'^) + m^(l + x) + 2mQm^{l - x) 



2{Q-^{l-x) + ml) 

where the splitting kernel p'^J reads 

pfj{x) = l-2x + 2x^. 



(A.16) 



(A.17) 



+ C>(e) (A.18) 



(A.19) 



B. Phase space mappings 

In this appendix we give the various phase space mappings used in the construction of our 
subtraction terms. These mappings must (a) conserve four momentum and (b) maintain 
the on-shellness of the particles involved. Three configurations need to be considered: final- 
final, initial-final and initial-initial, and depending if the mapping is related to a single or 
a double unresolved configuration 3 — > 2 or 4 — > 2 mapping procedures are needed. 

For processes involving massive particles the mapping procedures remain unaltered. 
As mentioned in section 2, provided all variables and momenta are expressed in terms of 
invariants of the form Sij = 2pi ■ pj, massless phase mappings can be used for processes 
with massive particles. For final- final configurations we use the mappings given in [68], 
whereas initial- final and initial- initial mappings are taken from [79]. 
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B.l Single unresolved mappings 

Subtraction terms built with three-parton antennae, namely those present at NLO and 
those related to the do-^^^^ subtraction terms, need the definition of {3 2} phase space 
mappings in all three configurations. 

B.1.1 Final-final configurations 

In this case, the mapping can be symbolically understood as {i,j,k) {ijijk), and it is 
defined by 

pf^- = xp\ + rpj + zp)^ 

P^, = {1- x)pf + (1 - r)p^ + (1 - (B.l) 



where 



X = 



[(1 + p)Sijk - 2rsjk] 



p2 4r(l - r)sijSjk 



(B.2) 

This mapping preserves momentum conservation 

Pij + Pjk =Pi+ Pj + Pk (B.3) 
and satisfies the following properties 

Pij Pi Pjk Pk when j 0, 

Pij -^Pi+ Pj Pjk Pk when i\ \j, 

Pij Pi Pjk Pj + Pk when j\\k, (B.4) 

which guarantee the correct subtraction of infrared singularities. 

B.l. 2 Initial-final configurations 

For these cases the mappings are of the form k) — t- (i,jk), where the bar on the initial 
state momentum indicates that it is only rescaled. The remapped momenta are given in 
terms of the original ones by 

-u u 
Pi = XiP, 

P^k=Pj+Pk-i^-Xi)pt^ (B-5) 

Sij + Sjk — Sjk 
Sij + Sik 



where 



-59- 



This mapping satisfies the following properties, which are needed for the correct subtraction 
of infrared singularities: 



XiPi Pi Pjk Pk when j 0, 

XiPi Pi Pjk Pj + Pk when j\\k, 

XiPi ^ (1 - Zi)pi Pjk Pk when pj = ZiPi. (B.7) 

B.1.3 Initial-initial configurations 

The mapping used for these configurations is of the form k, ...,m, n, ...} — > {i, k, fh, h, ...}. 
It is explicitely given by 

-U u 
Pi = XiPi 



Pj = XkPk 



(B.8) 



where 

q'=pt+I^k-Pj^ r = pf+pt- (B.9) 

The particular feature about initial-initial mappings is that all final state momenta, even 
those that are not actually part of the antenna, need to be remapped in order to preserve 
momentum conservation. This is due to the fact that q lies along the beam axis while q is 
in general not along the beam axis. The rescaling variables Xi and are given by 



Xi — 



Sik Sjk / Sik Sij Sjk 



^ik ^ij V ^ik 



= J ^ Sij Sik Sjk Sij ^ ^^^^^ 

y Sik Sjk V Sik 

and the properties satisfied by this mapping which ensure the correct subtraction of infrared 
singularities are 

XiPi Pi XkPk -> Pk when j -)• 0, 

XiPi Zi)pi XkPk Pk when pj = ZiPi, 

XiPi Pi XkPk Zk)Pk when pj = ZkPk- (B.H) 



B.2 Double unresolved mappings 

The use of four-parton antenna functions in the da^'^^Q parts of the subtraction terms 
requires the definition of adequate {4 — t- 2} mappings which smoothly interpolate the 
kinematics of the different unresolved configurations. 
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B.2.1 Final-final configurations 

In this case the mapping is {i,j,k,l) {ijk,jkl), and the remapped momenta are given 

by 



with 



n = 

r2 = 
X = 



Pijk = ^Pi + ^iPj + ^"^Pk + ^Pl 

= (1 - + (1 - + (1 - + (1 - 



Sij + Sj}; + S^; 

Ski 

Sik + Sjk + Ski 



2(sij + Sik + Sii) 
+(ri - r2) 
1 

2(sii + Sji + Ski 
-{ri - r2 
(ri - r2^^ 



(1 + p)sijki - ri{sjk + 2sj7) - r2isjk + 2sm) 

Sij Ski SikSkl 
Sil 

(1 - - n{sjk + 2sij) - r2{sjk + 2sik) 

I SijSkl ~l" SikSjl 
Sil 



1 + 



„2 2 
^il^ijkl 



A(Sjj, S/j;, Sj;, Sjki S^ki Sjl) 



2 (ri(l - r2) + r2(l - n)) (sijSfci + SifeSji - SjkSn) 



+4ri(l - ri)sijSji + 4r2(l - r2)sikSki 

As usually, the function A('U, is defined as 

A('u, lu) = + + ■i/;^ — 2(n'u + uw + wit;). 
The behaviour of these mappings in the different double unresolved limits is 



Pijk Pi 
Pijk Pi+ Pj 
Pijk -^Pi+ Pj + Pk 

Pijk ^ Pi 
Pijk Pi 
Pijk -^Pi+ Pj 



Pjkl Pi 
Pjkl Pi 
Pjkl Pi 

Pjkl Pj +Pk+Pl 
Pjkl -^Pl+Pk 
Pjkl -^Pl+Pk 



(B.12) 



(B.13) 



(B.14) 



when J, — 7- 

when /c — )• and 

when i||j||/c 

when |/ 

when J —7- and k\\l 

when and k\\l. (B.15) 



This behaviour ensures that infrared singularities are properly subtracted. Furthermore, it 
can be shown that in single unresolved limits, the momentum mappings defined in cq.(B.12) 
collapse into an NLO 3—7-2 mapping. This allows the subtraction of single unresolved limits 
in four-parton antenna functions with products of three-parton antennae. 
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B.2.2 Initial-final configurations 

For initial-final subtraction terms the mappings are {i,j,k,l) {i,jkl), and the remapped 
momenta are given by 

-ri a 

= Xip7^ 

with the rescaling variable 

Sij "I" Sik ~^ Sil Sjk Sjl Ski ^ „\ 

Xi — — — . [^■^') 

Sij ~r Sik ~r Sil 

The properties satisfied by these mappings in double unresolved limits are 

XiPi Pi Pjki Pi when j,k^ 

XiPi Pi Pjki -^Pk+Pl when j -)• and k\\l 

XiPi (1 - Zi)pi Pjki pi when k ^0 and pj = ZiPi 

XiPi Pi Pjki ^ Pj + Pk + Pi whenjplll 

XiPi ^ {I- Zi)pi Pjki Pi when pj + Pk = ZiPi 

XiPi (1- Zi)pi Pjki Pk+Pl when pj = ZiPi and k\\l. 

(B.18) 

As before, it can be shown that in single unresolved limits eq.(B.16) collapses into an NLO 
mapping. 

B.2.3 Initial-initial configurations 

Like in the single unresolved case, phase space mappings for these present configurations 
consist of a rescaling of the initial state hard radiators momenta, and a boost of all final 
state momenta that preserves overall momentum conservation. Concretely the mappings 
are of the form k, I, ...,m, ...) — )• (i, T, m, ...), and the remapped momenta are 

-u u 

Pi = XiPi 

-a a 
Pi = XipI 

where = +Pj + + Pi and = p^ -\- pf. The rescahng variables are given by 



Xi 



Sil — Sjl — Ski 


1 Sij Sik Sil "1" Sjk Sjl Ski 


Sil - ^ij - ^ik V 


Sil 


Sil Sij Sik 


1 Sij Sik Sil "1" Sjk Sjl Ski 


Sil + Sjl + Ski V 


Sil 



y Sil + Sjl + Ski V Sil 

In the double unresolved limits the mappings in eq.(B.19) give 

XiPi Pi xipi pi when j,k^O 
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XiPi Pi 
XiPi (1 
XiPi (1 



xiPi (1 



zi)pi 



when J —5- and pk = zipi 
when k ^ Q and pj = ZiPi 
when Pj +pk = ZiPi 
when Pj +pk = zipi 
when = ZiPi and = zipi 



Zi)pi 
Zi)pi 



xiPi Pi 



xiPi Pi 
xiPi (1 
xiPi (1 



XiPi Pi 
XiPi (1 



Zi)pi 



zi)pi 
zi)pi 



(B.21) 



while in single unresolved limits they collapse into NLO mappings. 
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